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Study of the neutrinoless double beta decay and searches for the manifestation of the neutrino 
mass in ordinary beta decay are the main sources of information about the absolute neutrino mass 
scale, and the only practical source of information about the charge conjugation properties of the 
neutrinos. Thus, these studies have a unique role in the plans for better understanding of the whole 
fast expanding field of neutrino physics. 
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I. EXECUTIVE SUMMARY 



The physics addressed by this research program seeks to answer many of the Study's questions: 

1. Are neutrinos their own anti-particles? 

2. What are the masses of the neutrinos? 

3. Do neutrinos violate the symmetry CP? 

4. Are neutrinos the key to the understanding of the matter-antimatter asymmetry of the Universe? 

5. What do neutrinos have to tell us about the intriguing proposals for new models of physics? 

Only the research covered within this working group can answer the first and second of these fundamental questions. 
Among the ways to measure the neutrino mass, three are notable because they are especially sensitive: double-beta 
decay, tritium beta decay, and cosmology. Consequently, we have focused our report and recommendations on them. 

• Observation of the neutrinoless double-beta decay (Oz//?/3) would prove that the total lepton number is not 
conserved and would establish a non vanishing neutrino mass of Majorana nature. In other words, observation of the 
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Oi//3/? decay, independent of its rate, would show that neutrinos, unhke aU the other constituents of matter, are their 
own antiparticles. There is no other reahstic way to determine the nature - Dirac or Majorana, of massive neutrinos. 
This would be a discovery of major importance, with impact not only on this fundamental question, but also on 
the determination of the absolute neutrino mass scale, and on the pattern of neutrino masses, and possibly on the 
problem of CP violation in the lepton sector, associated with Majorana neutrinos. There is a consensus on this basic 
point which we translate into the recommendations how to proceed with experiments dedicated to the search of the 
QvPP decay, and how to fund them. 

• To reach our conclusion, we have to consider past achievements, the size of previous experiments, and the 
existing proposals. There is a considerable community of physicists worldwide as well as in the US interested in 
pursuing the search for the decay. Past experiments were of relatively modest size. Clearly, the scope of future 
experiments should be considerably larger, and will require advances in experimental techniques, larger collaborations 
and additional funding. In terms of the effective neutrino Majorana mass that can be extracted from the 
observed Qv^P decay rate, there are three ranges of increasing sensitivity, related to known neutrino-mass scales of 
neutrino oscillations. 

• The ^100-500 meV (m/3/3) range corresponds to the quasi- degenerate spectrum of neutrino masses. The motivation 

for reaching this scale has been strengthened by the recent claim of an observation of decay in ^^Ge; a claim 

that obviously requires further investigation. To reach this scale and perform reliable measurements, the size of the 
experiment should be approximately 200 kg of the decaying isotope, with a corresponding reduction of the background. 

This quasi-degenerate scale is achievable in the relatively near term, ^ 3-5 years. Several groups with considerable 
US participation have well established plans to build ~ 200-kg devices that could scale straight-forwardly to 1 ton 
(Majorana using ^^Ge, Cuore using ^^'^Te, and EXO using ^^^Xe). There are also other proposed experiments 
worldwide which offer to study a number of other isotopes and could reach similar sensitivity after further R&D. 
Several among them [e.g. Super-NEMO, MOON) have US participation. 

By making measurements in several nuclei the uncertainty arising from the nuclear matrix elements would be 
reduced. The development of different detection techniques, and measurements in several nuclei, is invaluable for 
establishing the existence (or lack thereof) of the decay at this effective neutrino mass range. 

• The ~20-55 meV range arises from the atmospheric neutrino oscillation results. Observation of {mpis) at this mass 
scale would imply the inverted neutrino mass hierarchy or the normal-hierarchy v mass spectrum very near the quasi- 
degenerate region. If either this or the quasi-degenerate spectrum is established, it would be invaluable not only for the 
understanding of the origin of neutrino mass, but also as input to the overall neutrino physics program (long baseline 
oscillations, search for CP violations, search for neutrino mass in tritium beta decay and astrophysics/cosmology, etc.) 

To study the 20-50 meV mass range will require about 1 ton of the isotope mass, a challenge of its own. Given the 
importance, and the points discussed above, more than one experiment of that size is desirable. 

• The ^2-5 mcV range arises from the solar neutrino oscillation results and will almost certainly lead to the Ov/Sf] 
decay, provided neutrinos are Majorana particles. To reach this goal will require ~100 tons of the decaying isotope, 
and no current technique provides such a leap in sensitivity. 

• The qualitative physics results that arisc^ from an observation of Oz//3/3 decay are profound. Hence, the program 
described above is vital and fundamentally important even if the resulting {mpf}) would be rather uncertain in value. 
However, by making measurements in several nuclei the uncertainty arising from the nuclear matrix elements would 
be reduced. 

• Unlike double-beta decay, beta-decay endpoint measurements search for a kinematic effect due to neutrino mass 
and thus are "direct searches" for neutrino mass. This technique, which is essentially free of theoretical assumptions 
about neutrino properties, is not just complementary. In fact, both types of measurements will be required to fully 
untangle the nature of the neutrino mass. Excitingly, a very large new beta spectrometer is being built in Germany. 
This KATRIN experiment has a design sensitivity approaching 200 meV. If the neutrino masses are quasi-degenerate, 
as would be the case if the recent double-beta decay claim proves true, KATRIN will see the effect. In this case 
the Oi//3/3-decay experiments can provide, in principle, unique information about CP-violation in the lepton sector, 
associated with Majorana neutrinos. 

• Cosmology can also provide crucial information on the sum of the neutrino masses. This topic is summarized in 
a different section of the report, but it should be mentioned here that the next generation of measurements hope to 
be able to observe a sum of neutrino masses as small as 40 meV. We would like to emphasize the complementarity of 
the three approaches, Oi//3/3 , (3 decay, and cosmology. 

Recommendations: 

We conclude that such a double-beta-decay program can be summarized as having three components and our 
recommendations can be summarized as follows: 
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1. A substantial number (preferably more than two) of 200-kg scale experiments (providing the capability to make 
a precision measurement at the quasi-degenerate mass scale) with large US participation should be supported 
as soon as possible. 

o Each such experiment will cost approximately $10M-$20M and take 3-5 years to implement. 

2. Conciu'rently. the development toward ^1-ton experiments (i.e. sensitive to ■\/Amf^) should be supported, 
primarily as expansions of the 200-kg experiments. The corresponding plans for the procurement of the enriched 
isotopes, as well as for the development of a suitable underground facility, should be carried out. The US funding 
agencies should set up in a timely manner a mechanism to review and compare the various proposals for such 
experiments which span research supported by the High Energy and Nuclear Physics offices of DOE as well as 
by NSF. 

o Each such experiment will cost approximately $50M-$100M and take 5-10 years to implement. 

3. A diverse R&D program developing additional techniques should be supported. 

The total cost of this described program will be approximately $250M over a 10 year period. 

• In addition to double-beta decay, other techniques for exploring the neutrino mass need to be pursued also. We 

summarize these recommendations as follows. 

1. Although the KATRIN is predominately a European effort, there is significant US participation. The design and 
construction of this experiment is proceeding well and the program should continue to be strongly supported. 

2. Research and development of other techniques for observing the neutrino mass kinematically should be encour- 
aged. 
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II. INTRODUCTION 

The standard model of electroweak interactions, developed in the late 1960's, incorporated neutrinos as left-handed 
massless partners of the charged leptons. The discovery of the third generation of quarks and leptons completed the 
model, and made it possible, in addition, to incorporate also a description of CP violation. Later efforts to unify 
the strong and electroweak interactions led to the development of Grand Unified Theories which provided a natural 
framework for neutrino masses, and motivated many experiments in the field. Studies of e+e^ annihilation at the 
Z-resonance peak have determined the invisible width of the Z boson, caused by its decay into unobservable channels. 
Interpreting this width as a measure of the number of neutrino active flavors, one can, quite confidently, conclude 
that there are just three active neutrinos with masses of less than Mz/"^- 

In parallel, the understanding of big-bang nucleosynthesis and the discovery of the cosmic microwave background 
illustrated the important role of neutrinos in the history of the early universe. Those developments also led to the 
possibility that neutrinos, with small but finite mass, could explain the existence of dark matter. Although it now 
appears that neutrinos are not a dominant source of dark matter in the universe, the experimental evidence obtained in 
the last decade for finite neutrino masses and mixing between generations is strong and compelling. Those discoveries, 
the first evidence for 'physics beyond the Standard Model' gives us an intriguing glimpse into the fundamental source 
of particle mass and the role of flavor in the scheme of particles and their interactions. The scale of neutrino mass 
differences motivates experimental searches for the neutrinoless double beta decay and end-point anomalies in beta 
decay described in this report. 

Despite the recent triumphs of neutrino physics, several fundamental questions remain to be answered to advance the 
field itself and its impact in general on the whole particle and nuclear physics, as well as astrophysics and cosmology. 
The studies of neutrinoless double beta decay and end-point anomalies in beta decay, in particular, are essential 
and unique in their potential to answer the first two of them and plays an important and equally unique role in the 
remaining ones: 

• Are neutrinos their own anti-particles? 

• What are the masses of the neutrinos? 

• Do neutrinos violate the symmetry CP? 

• Are neutrinos the key to the understanding of the matter-antimatter asymmetry of the Universe? 

• What do neutrinos have to tell us about the intriguing proposals for new models of physics? 

The present report is structured as follows: In this introductory section we provide the "Goal of the field" statement 
first in which we stress the fundamental importance of the distinction between the Dirac and Majorana neutrinos, 
and its relation to the existence of the neutrinoless double beta decay {Oiy(3f3). 

In the next section "Oi//3/3 and /3 decays and Oscillations" we briefly summarize the status of the neutrino oscillation 
studies and the values of the corresponding mass differences Arn^^ and mixing angles. Next we discuss the relations and 
constraints provided by the results of oscillation studies and the neutrino mass parameters extracted from 0^(3(3 and 
/3 decay experiments. We also show how this research fits into the larger picture of the whole neutrino field described 
in the other reports of this APS study. Moreover, we stress the importance of Majorana neutrinos for leptogenesis, 
i.e. for the explanation of the (tiny) baryon/photon ratio and the present excess of baryons over antibaryons. 

For Qv(}j3, the process that is observable only in heavy nuclei, the understanding of the nuclear structure plays 
an essential role in extracting the neutrino effective mass from the observed rate. We discuss the nuclear structure 
aspects in the section titled "Nuclear Structure Issues" and comment on the existing uncertainties as well as on the 
prospects of reducing them. 

The rest of the report deals with the experiments. In the section titled "Experimental Prospects for /3/3" we 
summarize the situation in experimental fijd and briefly describe the numerous proposals. A similar treatment for /3 
decay is discussed in "Experimental Prospects for (3 Decay" . The relationship between cosmology and neutrino mass 
is described in the final section "Cosmology and Neutrino Mass" . We end with a conclusion section. 

Throughout the report we often use results reported in earlier reviews by some of us Q, H, H, Q . 

At present, we do not know the absolute scale of the neutrino mass. There is an upper limit, of ^ few eV from 
combining the limits from the tritium (3 decay with the Am^ values from the oscillation studies. For some, but not 
all, of the neutrinos there is also a lower limit, simply vKm? . These limits show that neutrinos, while massive, are 
very much lighter than the other fundamental constituents of matter, the charged leptons and quarks. While we do 
not understand the mass values of any fermion, the huge difference in masses of neutrinos and all charged fermions 
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clearly requires an explanation. The usual one, like the see-saw mechanism, ties the neutrino mass with some very 
high mass scale. It also suggests that neutrinos, unlike all other fermions, are Majorana particles, i.e. they are their 
own antiparticles. 

The research discussed here, if successfull, would show whether these ideas are true or not. As stated already in 
the executive summary, observation of the Qv[3(3 decay would prove that the total lepton number is not conserved 
and would establish a nonvanishing neutrino mass of Majorana nature. In other words, observation of the Qv[3l3 
decay, independently of its rate, would show that neutrinos, unlike all the other constituents of matter, are their own 
antiparticles. There is no other realistic way to determine the nature - Dirac or Majorana, of massive neutrinos. This 
would be a discovery of major importance, comparable to the already discovered oscillations of atmospheric, solar 
and reactor neutrinos, and as important as a discovery of CP violation involving neutrinos. It would have impact 
not only on this fundamental question, but also on the determination of the absolute neutrino mass scale, and on 
the pattern of neutrino masses, and possibly on the problem of CP-violation in the lepton sector, associated with 
Majorana neutrinos. 

At the same time, beta-decay endpoint measurements search for a kinematic effect due to neutrino mass and 
therefore are frequently referred to as "direct searches" for neutrino mass. This technique, which is essentially free of 
theoretical assumptions about neutrino properties, is not just complementary to the search of Qv[3l3 decay. In fact, 
both types of measurements will be required to fully untangle the nature of the neutrino mass. 

The following sections describe the status of the field, and plans for further experiments. Determining the absolute 
neutrino mass scale, and finding whether neutrinos are indeed Majorana particles and thus that the lepton number is 
not conserved, would represent a major advance in our understanding of particle physics. 

Qvl3l3 AND P DECAY AND OSCILLATIONS 
Status of oscillation searches 

As is well known, the concept of neutrino oscillations is based on the assumption that the neutrinos of definite flavor 
(i/e, v^^Vt) are not necessarily states of a definite mass z^i, 1^3.... Instead, they are generally coherent superpositions 
of such states, 

i 

When the standard model is extended to include neutrino mass, the mixing matrix U is unitary. As a consequence 
the neutrino flavor is no longer a conserved quantity and for neutrinos propagating in vacuum the amplitude of the 
process — > vii is 

A{vi ^vt')^Yl Uue-'^U;,, , (2) 

i 

The probability of the flavor change for £ ^ £' is the square of this amplitude, -^1^1') = \A{v^ '^£')\'^- It is 

obvious that due to the unitarity of U there is no flavor change if all masses vanish or are exactly degenerate. The 
idea of oscillations was discussed early on by Pontecorvo 5, 6] and by Maki, Nakagawa and Sakata [7|. Hence, the 
mixing matrix U is often associated with these names and the notation Umns or Up]\{NS is used. 

The formula for the probability is particularly simple when only two neutrino flavors, i^i and Vf, mix appreciably, 
since only one mixing angle and two neutrino masses mi,mj are then relevant, 

P(f£ vi'^i) = sin^ 29 SIT? 

where the appropriate factors of Ti and c were included. Here ^ra^ji = rUj — mf is the mass squared difference. Due to 
the large difference in the two observed Am? values, this simple formula adequately describes most of the experiments 
as of now. 

In general, the mixing matrix of 3 neutrinos is parametrized by three angles, conventionally denoted as 9i2, ^13, O23, 
one CP violating phase S and two Majorana phases ai,a2 [MlM^I- Using c for the cosine and s for the sine, the 
mixing matrix U is parametrized as 



(3) 
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C12C13 S12C13 si3 \ / e'"i/^ vi 




C12C23 - Sl2S23Sl3e'^ S23Ci3e^^ e*"^/^ V2 ■ (4) 



-S12C23 - Ci2S23Sl3e 

S12S23 - Ci2C23Si3e** -C12S23 - Si2C23Si3e*'' C23Ci3e"' / \ J/3 

The three neutrino masses rrii should be added to the parameter set that describes the matrix (0J, representing 
therefore nine unknown parameters altogether. 

The evidence for oscillations of solar (y,,) and atmospheric {v^ and D^) neutrinos is compelling and generally 
accepted. 

Evidence for oscillations of the solar Vf. have been reported first by the pioneering Davis et al. (Homestake) 
experiment It has been confirmed and reinforced later by Kamiokande, SAGE, GALLEX/GNO and Super- 

Kamiokande experiments 0, 0| . 

The SNO solar neutrino experiment 0, 0| , in which it is possible to separately determine the flux of v^, neutrinos 
reaching the detector (through the charged current reactions) and the flux of all active neutrinos (through the neutral 
current reactions), made the conclusion that solar neutrinos oscillate, inescapable. 

Independently, KamLAND reactor antineutrino experiment has shown that neutrinos oscillate as well 0, ^3 ■ 
Moreover, the oscillation parameters extracted from that experiment agree perfectly with those from the solar u,, 
experiments. This agreement, expected by the CPT-invariance, shows that the formalism of oscillations, including 
the matter efltects, is well understood. Based on the combined analysis of these data the parameters Am|i (including 
its positive sign) and Oq ~ 9i2 have been determined with a remarkable accuracy. 

Oscillations of the atmospheric v^i {v^j) have been most clearly observed in the Super-Kamiokande experiment. In 
particular the observed zenith angle dependence of the multi-GeV and sub-GeV ^-like events 0, 0| represents a 
compelling evidence. (Indication for the atmospheric neutrino oscillations, based mostly on the /x/e ratio, existed 
for a long time.) As is well known, the SK atmospheric neutrino data is best described in terms of dominant two- 
neutrino — !■ Vr [D^ — + Vt) vacuum oscillations with maximal mixing. The analysis thus fixes the parameters 
|Am3i| ~ |Am§2l and Oatm = 023 (since cos 6113 1). 

Finally, the remaining angle, ^13 remains unknown, but is constrained from above by the reactor neutrino CHOOZ 
and Palo Verde experiments |20|, |21| . 

The Table n summarizes the present status of knowledge of the oscillation parameters (assuming three mass eigen- 
states, i.e. disregarding the possible existence of sterile neutrinos). 

Thus, two of the three angles, and the two mass square differences have been determined reasonably well. The 
unknown quantities, accessible in future oscillation experiments (and discussed elsewhere in these reports) are the 
angle ^13 and the sign of the Am|2 ^ Am§]^. If that sign is positive, the neutrino mass pattern is called a normal mass 
ordering (mi < 7712 < m^) and when it is negative it is called inverted mass ordering (m^ < mi < m2). The extreme 
mass orderings, mi < m2 ^ m3 and m3 <^ mi < m2, are called the normal and, respectively, inverted hierarchies. In 
addition, the phase d governing CP violation in the flavor oscillation experiments remains unknown, and a topic of 
considerable interest. Determination of the CP phase S is again extensively discussed elsewhere in this report. 

The remaining unknown quantities, the absolute neutrino mass scale, and the two Majorana phases ai and a2 are 
not accessible in oscillation experiments. Their determination is the ultimate goal of Qvj3(3 and /3 decay experiments. 



TABLE I: Neutrino oscillation parameters determined from various experiments (2004 status) 

Parameter Value zbla Comment 
Amii 8.2t° j X lO-'^ eV^ 

6*12 32.3°+2/ For 6I13 = 

lAmial 2.0t° j x IQ-^ eV^ 
sin^ 2023 > 0.94 For 613 = 

sin2 26li3 <0.11 For Am^t„ = 2 X 10-^ eV^ 
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Oscillations and direct neutrino mass measurements 



Direct neutrino mass measurements are based on the analysis ol the kinematics of charged particles (leptons, pions) 
emitted together with neutrinos (flavor states) in various weak decays. The most sensitive neutrino mass measurement 
to date, involving electron type neutrinos, is based on fitting the shape of the beta spectrum (see section VI below). 
In such measurements the quantity 

= ^^|C/e.Pm2^ (5) 

is determined or constrained, where the sum is over all mass eigenvalues rrii,. that are too close together to be resolved 
experimentally. 

A limit on m^^ implies an upper limit on the minimum value m^^^.^ of all TOi^. , independent of the mixing parameters 
Uei'- < '^I'ej the lightest neutrino cannot be heavier than m^^. This is, in a sense, an almost trivial 

statement. 

However, when the study of neutrino oscillations provides us with the values of all neutrino mass-squared differences 
Amfj (including their signs) and the mixing parameters |J7eiP, and the value of m^^ has been determined in a future 
experiment, then the individual neutrino mass squares can be determined: 

iTT-l, = "^i^e ^ X! l^eipAml (Am-^ = ml. - ml^) . (6) 

i 

On the other hand, if only the absolute values |Amfj | are known (but all of them), a limit on m^^ from beta decay 
may be used to define an upper limit on the maximum value rn^^^^ of m^. : 

In other words, knowing |Ato? | one can use a limit on m^^ to constrain the heaviest active neutrino. 



Oscillations, Qvj3l3 decay, and neutrino mass 



The neutrinoless double beta decay. 



{Z,A) -> (Z + 2,A) + er 



(8) 



violates lepton number conservation. It can be recognized by its electron sum energy spectrum. Since the nuclear 
masses are so much larger than the decay Q value, the nuclear recoil energy is negligible, and the electron sum energy 
of the Oi//3/3 is simply a peak at Tei + Te2 = Q smeared only by the detector resolution. 

The 0vP(3 decay involves a vertex changing two neutrons into two protons with the emission of two electrons and 
nothing else. One can visualize it by assuming that the process involves the exchange of various virtual particles, e.g. 
light or heavy Majorana neutrinos, right-handed weak interaction mediated by the Wn boson, SUSY particles, etc. 
No matter what the vertex is, the OvPP decay can proceed only when neutrinos are massive Majorana particles ^2^. 

Of primary interest is the process mediated by the exchange of light Majorana neutrinos interacting through the 
left-handed V — A weak currents. The decay rate is then. 



[T°/2(0+ 



0+)]- 



m: 



GT 



9a 



(9) 



where G^" is the accurately calculable phase space integral, {mpp) is the effective neutrino mass, and M^^, 

are the nuclear matrix elements. The problems associated with the evaluation of the nuclear matrix elements are 

discussed in Section IV, where also corrections to the nuclear structure dependent part of the decay rate are discussed. 

If the QvPP decay is observed, and the nuclear matrix elements are known, one can deduce the corresponding {mfjp) 
value, which in turn is related to the oscillation parameters by 



IE 



(10) 
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where the sum is only over hght neutrinos (m^^ < 10 MeV). The Majorana phases ai were defined earher in Eq.Q. 
If the neutrinos Vi are CP eigenstates, ai is either or tt. Due to the presence of these unknown phases, cancellation 
of terms in the sum in Ea. HlU|l is possible, and {mpp) could be smaller than any of the m^. even if all neutrinos Vi 
are Majorana particles. 

We can use the values in Tableland express the {mpp) in terms of the three unknown quantities: the mass scale, 
represented by the mass of the lightest neutrino m,,^.^, and the two Majorana phases. In doing so, it is useful to 
distinguish three mass patterns: normal hierarchy (NH), toi < m2 -C ms (i.e. ruv^^^ = mi), invert ed hiera rchy (IH), 
TO3 ^ mi < m2 (i.e. rn^^.^ = ma), and the quasi-degenerate spectrum (QD) where mu^.^ 3> v^jAml^ as well as 
m^„,„ > V|Am^i|. 

In the case of normal hierarchy, and assuming that m^^^.^ = mi can be neglected, one obtains 

(m^/,)^^- I^A^sin^ 012 cos2 013 + ^^^A^sin^ 0136-*"^! . (11) 

For 013 = and the parameter values listed in Table {mpp)^-^ — 2.6 ± 0.3 meV. On the other hand, if 
tan0^3 > sin^ 012 y^|Am|77Am|7 ^ 0.06 a complete cancellation might occur and {mpp)^^ might be vanishingly 
small. However, if m^^^^ > then {mpfj)^^ may vanish even for 0i3 = 0, see Fig. ^ 

In the case of the inverted hierarchy, and again assuming that m^^^.^ = can be neglected, one obtains 



(m^^)^^ ~ y^lAmiJcos^ 0i3^1 - sin^ 20i2 sin^ . (12) 

Thus, if 013 = and for the parameter values listed in Tabled {mpp)^^ ~ 14—51 meV, depending on the Majorana 
phases. 

Finally, for the quasi-degenerate spectrum 

(mpfs)^" ~ moKcos^ 0i2e"i + sin^ 0i2e*"^) cos^ 0i3 + sin^ 0i3|. (13) 

If 013 = and for the parameters listed in Table HJ (m^^)'^^ ~ (0.71 ± 0.29)mo. 

Detailed discussion of the relation between the (mBp) and the absolute neutrino mass scale can be found in numerous 
papers (see, e.g. some of the more recent references 123. 125I l2^ Esf 

In Fig^we show the plot of {mpp) versus rn,y^.^ using the oscillation parameters in Table HJ and allowing for the 
maximum value of 0i3 and one a variations of them. One can clearly see the three regions (NH, IH, QD). Thus, 
determination of the {mpp) value would allow, in general, to distinguish between these patterns, and to determine a 
range of m^^^.^ . One should keep in mind, however, that there are caveats to this statement for the situations where 
the corresponding bands merge (e.g. the IH and QD near m^^.^ ~ 0.05 eV). 

Despite this caveats, obviously, if one can experimentally establish that {mpp) > 50 meV, one can conclude that 
the QD pattern is the correct one, and one can read off an allowed range of m^^^.^ values from the figure. (The sign 
of Amg]^ ~ Am|2 will remain undetermined in that case, however.) 

On the other hand, if {mpp) 20-50 meV only an upper limit for the m^^.^ can be established, and the pattern 
is likely IH, even though exceptions exist. However, if (and that is unlikely in a foreseable future) the value of m^^^.^ 
can be determined independently, the pattern can be resolved. 

Finally, if one could determine that (m^^) < 10 meV but nonvanishing (which is again is unlikely in a foreseable 
future), one could conclude that the NH pattern is the correct one. 

Altogether, observation of the Qufifi decay, and accurate determination of the {mpp) value would not only establish 
that neutrinos are massive Majorana particles, but would contribute considerably to the determination of the absolute 
neutrino mass scale. Moreover, if the neutrino mass scale would be known from independent measurements, one could 
possibly obtain from the measured (m^^) also some information about the CP violating Majorana phases. 

Absolute neutrino mass scale 

As shown above, 0v(3j3 and (3 decays both depend on different combinations of the neutrino mass values and 
oscillation parameters. The Qv(3(3 decay rate is proportional to the square of a coherent sum of the Majorana neutrino 
masses because {mpf}) arises from exchange of a virtual neutrino. On the other hand, in beta decay one can determine 
an incoherent sum because a real neutrino is emitted. 
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FIG. 1: Effective Majorana mass (m^^) versus the minimum mass iriu^-^. The different mass patterns are indicated. The 
shaded region corresponds to the best values of oscillation parameters, and S13 — 0. The dashed lines indicate the expanded 
range corresponding to the la errors of the oscillation parameters and the maximum allowed 613. Note that the uppermost 
line is unchaged (within this scale) in that case. 



Quite different source of information is based on cosmological and astrophysical observations where the density of 
the primordial neutrino sea is determined or constrained and thus a parameter proportional to the sum of the neutrino 
masses is determined. 

Massive neutrinos would contribute to the cosmological matter density an amount, 

Q^h^ ^J:m^j92.5 eY , (14) 

where il^ is the neutrino mass density relative to the critical density and lOO/i is the Hubble constant in km/s/Mpc. 
From the requirement that the neutrinos left over from the Big Bang do not overdose the universe an upper limit, 
with a minimum assumptions (essentially just the requirement of stability), is obtained 

46 eV 

< ^ , (15) 

where N^, is the number of neutrino species with standard weak interactions |29| . 

More restrictive limits are obtained from the requirement that excessive free streaming in the early universe would 
not suppress small scale power of the observed matter distribution. The relation between the damping scale dps 
caused by free streaming, and the neutrino mass is approximately 

dps (Gpc) ^ 1/m, (eV) . (16) 

The data on Cosmic Microwave Background (CMB) and large scale galaxy surveys can be used to constrain N^rrii, 
for the quasi-degenerate neutrino mass spectrum, and thus also for various assumed number of neutrino flavors 
N^. The following Table ITU is based on js^l- Different analyses with different assumptions typically reach similar 
conclusions, suggesting that these limits are fairly robust (see more discussion further in this report). 

For completeness, note that, in principle, neutrino mass can be also extracted from the time of flight determination 
of neutrinos from future galactic supernova. However, one does not expect to be able to reach sub-eV sensitivity with 
this method (see e.g.js^)- 

It is worthwhile to stress that the various methods that depend on the neutrino absolute mass scale are comple- 
mentary. If, ideally, a positive measurement is reached in all of them (0^/3/3 (3 decay, cosmology) one can test the 
results for consistency and perhaps determine the Majorana phases. We illustrate the idea |3 in Fig. |2 using a 
two- neutrino-species example of such a set of measurements. (A 3-species example is discussed in Ref. '3|.) We took 
the mixing matrix and Am^ to be the best fit to the solar-neutrino data, with an arbitrary value for the Majorana 
phase a (of which there is only one) of 2.5 radians. We then made up values for S, {mpp), and (m/3) assuming them to 
be the results of pretend measurements. Each curve in the vs. mi graph is defined by one of these measurements. 
We chose the value of S (from cosmology) to be 600 meV, corresponding to a quasidegenerate hierarchy, and let (to/j) 
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TABLE II: Limits on NvUii, and on assuming quasi-degenerate mass spectrum. In the last column is the heaviest neutrino 
mass assuming that one neutrino dominates the sum. 

iV,. N^^m^ (eV) (eV) m^„„^ (eV) 



3 


1.01 


0.34 


0.73 


4 


1.38 


0.35 


1.05 


5 


2.12 


0.42 


2.47 


6 


2.69 


0.45 


4.13 



= 300 meV and {mpp) — 171 meV. The TO2 versus mi curves from the "measurements" of the oscillation parameters, 
S, and P decay are: 

m2 = S — TOi 

m2 = \J m\ + 5m\2 

m2 = ^J{{'mp)lU,2?~{niiUeilU,2? ■ (17) 

The 13(3 constraint is a little more complicated than that from (3 decay. The curve is also an ellipse but rotated 
with respect to the axes. All of the equations express m2 in terms of mi and measured parameters and all should 
intersect at one (mi,m2) point. However, because the point is overdetermined, the (3[3 ellipse will intersect only for a 
correct choice of a. This provides a way to determine a. In Figure |2] we drew the /?/? ellipse for a = 2.0 radians and 
for the "true" value of 2.5 radians to show how the intersection does indeed depend on a correct choice of the phase. 

Many authors have examined the potential to combine measurements from (3(3 decay, tritium /3 decay, and cosmology 
to determine the Majorana phases. (See e.g. [stUssHs^ .I 



Leptogenesis 

Leptogenesis is a natural possibility when the see-saw mechanism is realized by Nature. Since this mechanism 
predicts the light neutrinos to be of Majoarana type, Qvl3(3 is also expected. Though the details of the connection 
of the high and low energy parameters remain to be investigated in a given model, leptogenesis always comes with 
QvPP. The only exception is when in the basis in which the charged leptons are diagonal the neutrino mass matrix 
has a zero entry in the ee element, resulting in no ncutrinoless double beta decay. This implied texture will however 
have interesting physical implications by itself. In this section we discuss the relationship between neutrino mass and 
leptogenesis. 

The origin of the matter-antimatter asymmetry is one of the most important questions in cosmology. The presently 
observed baryon asymmetry is 

= 6.5 X 10-1°. (18) 

s 

In 1967 A. Sakharov suggested that the baryon density can be explained in terms of microphysical laws Three 
conditions need to be fuUfilled: 

• Baryon number (or Lepton number, for the leptogenesis mechanism) violation. Assuming that the initial baryon 
asymmetry is negligible, as implied by inflationary models, baryon number violation is required to produce a 
net final baryon asymmetry. 

• C and CP violation. The CP symmetry is the product of charge conjugation and parity. If CP is conserved, 
every reaction which produces a particle will be accompanied by a reaction which gives its antiparticle, with no 
creation of a net baryon number. 

• Departure from thermal equilibrium. 

Several mechanism have been proposed to understand the baryon asymmetry. 
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FIG. 2: A consistency plot for the neutrino mass eigenvalues mi and m2, for various hypothetical measurements. This set of 
curves indicates how measured values of S, {mpp), Am^^i, and {mfj) constrain the mass eigenvalues. The Qvf3(3 curve has been 
drawn for an incorrect value of the phase in the bottom panel to indicate the sensitivity of this technique for extracting the 
CP-violating phase. See text for further description of the parameters used to draw the curves. 



- Planck-scale baryogenesis. While we can expect a quantum theory of gravity not to preserve any quantum number, 
these effects would be relevant only at very early times and the baryon asymmetry generated in this way would be 
subsequently diluted in the inflationary epoch. 

- GUT baryogenesis (see, e.g., ref. |2J|)- GUT theories assume the existence of an underlying theory at high 
scale with a simple gauge group which is then broken to the Standard Model gauge group at low energy scale. The 
unification scale is expected to be at ~ 10^^ GeV. In general these models have lepton and baryon number violation 
and have many new parameters with respect to the Standard Model which can provide new sources of CP-violation. 
However cosmological considerations, mainly overproduction of gravitinos, suggest that the reheating temperature 
after inflation should be not higher than ~ 10^ GeV. If so, any baryon asymmetry produced at the GUT scale would 
be later diluted and negligible. 

- Electroweak baryogenesis. The interactions of the Standard Model, even if preserving the lepton and baryon 
number at the perturbative level, do violate it non-perturbatively [HI. The ground states of a non-abelian gauge 
group are labeled by an integer n and are separated by energy barriers one from the other. At zero temperature, 
the instantons transitions between different vacua are exponentially suppressed, because they are due to quantum- 
tunneling effects. At high temperatures above the electroweak breaking scale, these transitions can proceed through 
sphaleron processes. This produces a net baryon number. It was pointed out the Standard Model can satisfy all 
the three Sakharov condition for baryogenesis |36|. The departure from equilibrium is realized at the electroweak 
symmetry breaking scale when the Higgs field acquires a vacuum expectation value different from zero and a phase 
transition takes place. However, the generated baryon asymmetry is by far too small to account for the observations. 
This is mainly due to the fact that the departure from equilibrium is not strong enough, as implied by a heavy Higgs, 
iTih > 80 GeV. 
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Supersymmetric extensions of the Standard Model contain new sources of CP-violation and an enlarged set of pa- 
rameters which might produce a first order phase transition. The parameter space which allows for a sizable baryon 
number generation is rather small (it requires in particular a light stop) and it will be soon tested. 

- Affleck-Dine baryogenesis |37| . Scalar fields carrying baryon and lepton number are naturally present in super- 
symmetric extension of the Standard Model (e.g. squarks and sleptons). A coherent scalar field, a field which has 
large vacuum expectation value, can in principle carry a large baryon number. This condensate can be produced 
in the Early Universe by the evolution of the scalar field along fiat directions of its potential. The decay of the 
condensate into ordinary matter is then responsible for the baryon asymmetry of the Universe, as the baryon charge 
of the condensate is transferred to quarks. This is still a viable explanation for the baryon asymmetry of the Universe. 

- Leptogenesis [l^- Let us notice that B — L is conserved both at the perturbative and non-perturbative level. 
This implies that if one creates a net B — L, (e.g., a lepton number), the sphaleron processes would leave both baryon 
and lepton number comparable to the original B — L. This idea is implemented in the leptogenesis scenario [Sg- 
Leptogenesis is particularly appealing because it takes place in the context of see-saw models |39l|, which naturally 
explain the smallness of neutrino masses. The see-saw mechanism requires the existence of heavy right-handed (RH) 
Majorana neutrinos, completely neutral under the Standard Theory gauge symmetry group. Consequently, they can 
acquire Majorana masses that are not related to the electroweak symmetry breaking mechanism, and can, in principle, 
be much heavier than any of the known particles. Introducing a Dirac neutrino mass term and a Majorana mass term 
for the right-handed neutrinos via the Lagrangian 



-C = 17n {mD)^J Nrj -t- - (TVffi)^ {Mr),j Nrj , (19) 

leads, for sufficiently large Mr, to the well know see-saw [39] formula 

ruv ^ —ruD M^^ , (20) 

= C^PMNS Dm t^PMNS (21) 

where terms of order 0{M^'^) are neglected. Here Dm is a diagonal matrix containing the masses toi,2,3 of the 
three light massive Majorana neutrinos and C/pMNS is the unitary Pontecorvo-Maki-Nagakawa-Sakata lepton mixing 
matrix. 

The CP-violating and out-of-cquilibrium decays of RH neutrinos produce a lepton asymmetry |38l | that can be 
converted into a baryon asymmetry through anomalous electroweak processes j3a,|40j. The requisite CP violating 
decay asymmetry is caused by the interference of the tree level contribution and the one-loop corrections in the decay 
rate of the three heavy Majorana neutrinos, Ni <f>^ and Ni 1^: 

_ T{N, £+) - T{Ni $+ j-) 

~ r(iV, -> i+)+T{N, ) 

1 , ^ Im(mt,mr))?. (2^) 



where $ and £ indicate the Higgs field and the charged leptons, respectively. Here w ~ 174 GeV is the electroweak 
syinmetry breaking scale and Xj = M'j/Mf. The functions / and g stem from vertex 0| and from self-energy 
1421 contributions. 



f{x)^^[l-{l + x) In(ii^)) 
^(-) - S-x 



(23) 



For X ^ 1, i.e., for hierarchical heavy Majorana neutrinos, one has f (x) + g{x) ~ Then, the baryon asymmetry 

is obtained via Yb = a{K/ g*) ei, where a ~ —1/2 is the fraction of the lepton asymmetry converted into a baryon 
asymmetry [s^ 9* — 100 is the number of massless degrees of freedom at the time of the decay, and k is an 
efficiency factor that is obtained by solving the Boltzmann equations Typically, one gets Yb ~ 6 x 10^^" when 
£i ~ (10~^ — 10""^) and k, ^ (lO^'^ — 10^^). A similar estimate of Yb is obtained in the supersymmetric (SUSY) 
theories as well [i^ . 

In non-supersyymetric models which embedd the see-saw mechanism, lepton flavour violating (LEV) charged lepton 
decays such as/z^e-l-7, r^/i-|-7, r— s-e-l-7, are predicted to take place with branching ratios which are strongly 
suppressed SUSY theories have additional sources of lepton charge non-conservation. If SUSY is broken above 
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the RH Majorana mass scale, as, e.g, in gravity-mediated breaking scenarios, there are renormalization group effects 
that gene rate new lepton charge non-conserving coupHngs at low energy even if such couplings are absent at the GUT 
scale |43,|4a|. In contrast to the non-supersymmetric case, the LVF processes can proceed with rates and cross-sections 
which are within the sensitivity of presently operating and proposed experiments. Under the assumption of flavour 
universality of the SUSY breaking sector (scalar masses and trilinear couplings) at the GUT scale Mx, using the 
leading-log approximation, one can estimate 0,0,0] the branching ratio of the charged lepton decays ii ^ ij +j, 
^T,fi,e for = 3, 2, 1, i > j, 



BR{1, I, + 7) 



{mu diag(log^,^,^) ml\, 



(3 + ao)mQ 

Jtt^ m% GpfJ (24) 



tan^ P , 



where mg denotes a slepton mass , m o is the universal mass scale of the sparticles and oq is a trilinear coupling (all 
at Mx) (for details, see, e.g., ref. |43|). 

Establishing a connection between the parameters at low energy (neutrino masses, mixing angles and CP-violating 
phases), measurable in principle in present and future experiments, and at high energy (relevant in leptogenesis) has 
gathered a great interest in the last few years. The number of parameters in the full Lagrangian of models which 
implement the see-saw mechanism is larger than the ones in the low-energy sector: in the case of 3 light neutrinos 
and three heavy ones, at high energy the theory contains in the neutrino sector 18 parameters of which 12 real ones 
and 6 phases, while at low energy only 9 are accessible - 3 angles, 3 masses and 3 phases. The decoupling of the 
heavy right-handed neutrinos implies the loss of information on 9 parameters. This implies that reconstructing the 
high energy parameters entering in the see-saw models from the measurement of the masses, angles and CP-violating 
phases of rrii, is in general difficult, if not impossible, and depends on the specific model considered. 

Using the weak basis in which both Mr and the charged lepton mass matrix are real and diagonal, we find useful 
to parametrize the Dirac mass by the biunitary or the orthogonal parametrizations. 

Biunitary parametrization. We can write the complex 3x3 Dirac mass as (see, e.g., ref 0|): 

mo = Ul m^'"s Ur , (25) 

where Ul and Ur are unitary 3x3 matrices and m^'^® is a real diagonal matrix. All the CP-violating phases are 
contained vhUl and Ur. 

Orthogonal parametrization. 

By using the see-saw formula, eq. (jSOJ, we can express m^i as [i^li^: 

mo - t C/PMNS Dll"" R M'J^ , (26) 

where Dm is the diagonal real matrix which containes the low-energy light neutrino masses, and i? is a complex 
orthogonal matrix. R contains 3 real parameters and 3 phases. 

The use of the two indicated parametrization clarifies the dependence of the processes we are interested in, e.g. 
leptogenesis and LFV charged lepton decays, on the different parameters entering in m^). In particular we have that: 
- for leptogenesis, the decay asymmetry ei depends on the hermitian matrix m\ymu: 

f C^fl(™D"')'f^fl - bi-unitary; 
m'jjTnD = < (27) 
[ M^^^ R^ D^RMr^ , orthogonal. 

We can notice that the PMNS unitary mixing matrix U does not enter esplicitly into the expression for the lepton 
asymmetry. 

- Concerning LFV decays, in supersymmetric versions of the see-saw mechanism, their rates depend approximately 
on H^iSj 

ul {mt')' Ul , 

mom^ = <^ (28) 



C/pMNS R Mr i?t dU^ C// 



PMNS ' 



by using the biunitary (upper expression) or the orthogonal parametrization (lower expression) . 
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- finally, in the biunitary parametrization, tiie neutrino mass matrix can be written as: 

TO. = -Ul m^'"s jj^ jjT ^diag _ ^29) 

This shows that the phases in JJpMNS receive contribution from CP-violation both in the right-handed sector, 
responsible for leptogenesis, and in the left-handed one, which enters in lepton flavour violating processes. Due to the 
complicated way in which the high energy phases and real parameters enter in m^,, eq. H29|) . if there is CP-violation at 
high energy, as required by the leptogenesis mechanism, we can expect in general to have CP-violation at low-energy, 
as a complete cancellation would require some fine-tuning or special forms of TOd and M^. Let us mention that an 
observation of CP-violation at low energy, however, does not imply necessarily CP-violation in U r as it might receive 
all its contributions from U l ■ 

More specifically, from eq. H29(l . we see that, in general, there is no one-to-one link between low energy CP-violation 
in the lepton sector and the baryon asymmetry: a measurement of the low energy CP-violating phases does not allow 
to reconstruct the leptogenesis phases [soj . However most specific models allow for such a connection. In particular 
if the number of parameter is reduced in tod, then a one-to-one correspondence between high energy and low energy 
parameters might be established. This can be achieved in models which allow for CP-violation only in the right- 
handed sector, that is in Ur, or which reduce the number of independent parameters at high energy, for example by 
requiring only two right-handed neutrinos |5lj . Each model of neutrino mass generation should be studied in detail 
separately to establish the feasibility of the leptogenesis mechanism [s^ . 

Furthermore the amount of baryon asymmetry depends on the type of li ght and heavy neutrino mass spectra and 
succesful leptogenesis imposes some constrains on the allowed mass scales |48[ 0] ■ For example in ref. the 
case of strongly hierarchical see-saw models was analyzed, assuming a hierarchical structure of the heavy Majorana 
neutrino masses and of the neutrino Dirac mass matrix mo- In order to produce a sufficient amount of baryon 
asymmetry via the leptogenesis mechanism, the scale of tod should be given by the up-quark masses. In this case the 
CP violation effects in neutrino oscillations can be observable, but, in general, there is no direct connection between 
the latter and the CP violation in leptogenesis. In ref. [3| it is also shown that, in the case of thermal leptogenesis 
with a hierarchical structure of the heavy Majorana masses, a strong bound on the light neutrino masses can be put 
of order 0.1 eV. This bound can be weakened if one assumes non-thermal right-handed neutrino production and/or 
an enhancement of the decay asymmetry due to quasi-degenerate right-handed neutrinos (resonant leptogenesis). For 
example, in ref. the case of quasi-degenerate spectrum of light neutrinos, toi^2,3 — "m^ > 0.2 eV, was considered. 
It was shown that leptogenesis is feasible and can provide the observationally required baryon asymmetry of the 
Universe. Furthermore it was pointed out that the branching ratios for the lepton flavour processes, n ^ e + 
T ^ II + ^ and r — !■ e -|- 7 decay, might be close to the present bounds and that they depend on the Majorana 
CP-violating phases which enter in (to) . 

The possible observation of Ou/SP decay would play an important role in understanding the origin of the baryon 
asymmetry as it would imply that lepton number (one of the main conditions for leptogenesis) indeed is not conserved. 
Furthermore the Majorana nature of neutrinos would be established: the see-saw mechanism would be regarded as 
a reasonable explanation of neutrino mass generation. Leptogenesis naturally takes place in this scenario. Finally 
the observation of CP-violation in the lepton sector, in neutrino oscillation experiments and/or Oi^/3/3 decay, would 
suggest the existence of CP-violation at high energy, which might be related to the one responsible for leptogenesis. 

In conclusion, the observation of lepton number violation in Qvf3j3 decay and, in addition, possibly of CP-violation 
in the lepton sector, would be a strong indication, even if not a proof (as it is not possible to reconstruct in a model 
indipendent way the high energy parameters from to.), of leptogenesis as the explanation for the observed baryon 
asymmetry of the Universe. 

Alternative tests of lepton number conservation 

Our main emphasis is on the study of the 0p(3P decay, 

{A,Z) — > {A,Z + 2) + e- + e- . (30) 

However, alternatively, similar physics could be tested in other processes, e.g. in 

+ {A,Z) e+ + {A,Z-2) (31) 

K+ — y fj.+ + ii+ + TT- (32) 
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which can be mediated by the exchange of virtual massive Majorana neutrinos |55|, |5g and where rather restrictive 
Umits on the corresponding branching ratios exist. These Lepton number violating decay rates are typically pro- 
portional to an average Majorana neutrino mass of the form: (m«')^ = | UiiUi'ie°'^mi\'^ , where ai/2 denotes the 
Majorana phases (which can lead to destructive interference). 

Searches for muon-positron conversion and rare kaon decays yield: (m^e) < 17(82) MeV, and {m^^j) < 

4 X 10^ MeV, respectively (see and references therein). These bounds therefore do not constitute a useful constraint 
for the mass pattern of the three light neutrinos. Similarly to the kinematic tests, new Majorana mass experiments 
are focusing on process 0v(3l3 offering by far the best sensitivity. 

Note that other lepton number violating processes might also exist, causing for example the emission of the De 
antineutrinos by the Sun. Even though the existence of such processes would again signal that neutrinos are massive 
Majorana particles, the interpretation of the rate in terms of neutrino masses is not straigthforward. Instead, other 
mechanisms might be at play, like the spin-flavor precession combined with neutrino oscillations (relevant parameter 
is then fi^Bx, where /ij^ is the transition magnetic moment and Bx is the magnetic field in the Sun |60j ) or neutrino 
decay (the relevant parameter is then r/m of the heavier neutrino oT]). 



Mechanism of Oi//3/3 decay 

Although the occurrence of 0^/3/3 decay implies the existence of massive Majorana neutrinos their exchange 
need not be the only and even dominant contribution to the decay rate. Almost any physics that violates the total 
lepton number can cause Qu(3(3 decay. A heavy Majorana neutrino can be exchanged, or supersymmetric particles, 
or a leptoquark. Right-handed weak currents, either leptonic or hadronic, can cause the absorption of an emitted 
virtual neutrino without the helicity flip that depends on the neutrino mass. Moreover, heavy particle exchange, 
heavy neutrinos, supersymmetric particles (see, e.g., 1^3 )j resulting from lepton number violation dynamics at some 
scale above the electroweak one can lead to Qvj3(3 decay with the same single electron spectra as in the favored case of 
the light Majorana neutrino exchange. These alternative mechanisms, if dominant, would not allow (or would make 
it very difficult) to extract the effective Majorana mass {mpp) from the measured Ovfifi decay rate. 

The relative size of heavy (Ah) versus light particle [Al) exchange contributions to the decay amplitude can be 
crudely estimated as follows |6.ll |: 

- Gf^. Ah - G^ — , — ^ , (33) 

where k ~ (50 MeV) is the typical light neutrino virtual momentum, and A is the heavy mass scale relevant to the 
lepton number violation dynamics. Therefore, Ah/Al ~ 0(1) for to/j^ ^ 0.1 — 0.5 eV and A '--^ 1 TeV, and thus 
the TeV scale leads to similar Oj^/3/3 decay rate as the exchange of light Majorana neutrinos with the effective mass 
TO/3/3 ^ 0.1 — 0.5 eV. (Note that if this heavy mass scale A is larger than 3-5 TeV the Al will clearly dominate.) 

The exchange of heavy particles involves short-range propagators that give rise to decay rates of the same form as in 
the mass-driven mode: simple two-s- wave-electron phase space multiplied by the square of an amplitude. Two-nucleon 
correlations do not suppress the effects of heavy particles, which can be transmitted between nucleons by pions |62j . 
The angular distributions and single-electrons spectra will therefore be the same in all these processes. One way to 
distinguish one from another is to take advantage of the different nuclear matrix elements that enter the amplitudes 
(leading to different total decay rates). Unknown parameters such as the effective neutrino mass or the trilinear 
i?-parity-violating supersymmetric coupling (violation of R parity naturally accompanies Majorana neutrino-mass 
terms) also enter the rates, so several transitions would have to be measured. This might be best accomplished [H^ 
by measuring transitions to several final states in the same nucleus, but if the matrix elements can be calculated 
accurately enough one could also measure the rates to the ground states of several different nuclei. The problems in 
determining the source of Qv(3f3 decay are mitigated by constraints from other experiments on many extra-Standard 
models. Some of these constraints will be much stronger once the Large Hadron Collider comes on line. 

Alternatively, it has been argued in [g^ that the study of the lepton flavor violating processes ^ —> e conversion, 
and ^ — + e 4- 7 will provide important insight about the mechanism of the OvP/] decay. In particular, the ratio of 
the branching fractions of these processes provides a diagnostic tool for establishing whether the exchange of light 
Majorana neutrinos is the dominant mechanism (and hence determination of (TO/3/3) is possible) or whether further 
more involved analysis is needed. 
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III. NUCLEAR STRUCTURE ISSUES 

The observation of Qv(5[3 decay would immediately tell us that neutrinos are Majorana particles and give us an 
estimate of their overall mass scale. But without accurate calculations of the nuclear matrix elements that determine 
the decay rate it will be difficult to reach quantitative conclusions about masses and hierarchies. 

Theorists have tried hard to develop many-body techniques that will allow such calculations. In order to test the 
calculations they have tried to calibrate them against related observables: 2v(3(i decay, ordinary /3+ and f3~ decay, 
Gamow- Teller strength distributions, odd-even mass differences and single-particle spectra. They have tried to exploit 
approximate isospin and 5L/(4) symmetries in the nuclear Hamiltonian and to extend well-known many-body methods 
in novel ways. In spite of all this effort, we know the matrix elements with only limited accuracy. In this section we 
review the state of the nuclear-structure calculations and discuss ways to improve them. 

Most recent attempts to calculate the nuclear matrix elements have been based on the neutron-proton Quasiparticle 
Random Phase Approximation (QRPA) or extensions to it. Of those that haven't, the most prominent are based on 
the shell model (SM). While the two methods have much in common — their starting point is a Slater determinant of 
independent nucleons — the kinds of correlations they include are complementary. The QRPA treats a large fraction 
of the nucleons as "active" and allows these nucleons a large single-particle space to move in. But RPA correlations 
are of a specific and simple type best suited for collective motion. The shell model, by contrast, treats a small fraction 
of the nucleons in a limited single-particle space, but allows the nucleons there to correlate in arbitrary ways. That 
these very different approaches yield similar results indicates that both capture most of the important physics. That, 
by itself, is encouraging and restricts the possible values of nuclear matrix elements considerably. 

QRPA 

The application of QRPA to (3(3 decay began with the realization by [g^ that in the QRPA the neutron-proton 
{np) particle-particle (i.e. pairing-like) interaction, which has little effect on the collective Gamow- Teller resonance, 
suppresses 2v(3(3 rates considerably. Soon afterward, ^ and [z3 

demonstrated a similar though smaller effect on 
Qv(3(3 decay. It was quickly realized, however, that the QRPA was not designed to handle realistic np pairing; the 
calculated half-lives were unnaturally sensitive to the strength of the pairing interaction. As a result, the rates of 
(3(3 decay, particularly 2v(3(3 decay, were hard to predict precisely because a small change in a phenomenological 
parameter (the strength of np isoscalar pairing) caused a large change in the lifetimes and eventually the breakdown 
(called a "collapse" ) of the entire method when the parameter exceeds some critical value. Most recent work in the 
QRPA has aimed at modifying the undesirable aspects of the method so that its sensitivity to np pairing becomes 
more realistic. 

There has been a large number of attempts to extend the validity of QRPA. We will not list them here, or discuss 
them in detail. Gomprehensive discussion can be found in Q . Earlier and even more complete discussion of the issues 
involved can be found in |7ll |. 

Unfortunately, it's not clear which of the extensions of the standard QRPA is best. Some of them violate an 
important sum rule for single (3 strength, and studies in solvable models suggest that the reduced dependence (at 
least of the 'renormailzed QRPA (RQRPA)) on neutron-proton pairing may be spurious resulting from an artificial 
reduction of isoscalar pairing correlations. There are not many obvious reasons to prefer one of these QRPA extensions 
or the other. 

Recently, one of the methods, RQRPA, was modified even further, so that the single (3 sum rule was restored. 
The resulting method, called the "Fully Renormalized QRPA" has yet to be applied to Qv(3(3 decay. Even more 
recently, Ref. [t^ raised the issue of nuclear deformation, which has usually been ignored in QRPA-like treatments 
of nearly spherical nuclei. (Psuedo-SU(3)-based truncations have been used to treat it in well-deformed nuclei ^3|)- 
The authors of [72,J argued that differences in deformation between the initial and final nuclei can have large effects 
on the 2v(3(3 half-life. These ideas, too, have not yet been applied to Qv(3(3 decay. 

The profusion of RPA-based extensions is both good and bad. The sheer number of methods applied gives us a kind 
of statistical sample of calculations, which could give an idea of the theoretical uncertainty in the matrix elements. 
But the sample may be biased by the omission of non-RPA correlations in all but a few calculations. Also, different 
calculations are done with different ways of fixing the parameters, e.g. some insist to reproduce the 21/(3(3 rate, other 
are unable to do that (or do not insist on it), etc. It is not clear, therefore, that the spread of the calculated values 
indeed corresponds to the spread of the possible matrix elements within even just QRPA and its extensions. 
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Shell Model 



The obvious alternative to QRPA, and the current method of choice for nuclear structure calculations in heavy 
nuclei where applicable, is the shell model. It has ability to represent the nuclear wave function to arbitrary accuracy, 
provided a large enough model space is used. This caveat is a huge one, however. Current computers allow very large 
bases (millions of states), but in heavy nuclei this is still not nearly enough. Techniques for constructing "effective" 
interactions and operators that give exact results in truncated model spaces exist but are hard to implement. Even 
in its crude form with relatively small model spaces and bare operators, however, the shell model seems to offer some 
advantages over the QRPA. Its complicated valence-shell correlations, which the QRPA omits (though it tries to 
compensate for them by renormalizating parameters) apparently affect the 13(3 matrix elements |74|. 

The first modern shell-model calculations of (3(3 decay date from the work Haxton and Stephenson vS] and references 
therein. Only a few truly large-scale shell model calculations have been performed. The heavy deformed (3(3 nuclei, 
"^^^U , and ^^"Nd, for example, require bases that are too large to expect real accuracy. Realistic work has thus been 
restricted to few nuclei, in particular to ''^Ca, "^^Ge, and ^^^Xe, though less comprehensive calculations have been 
carried out in several other nuclei [7ll | . 

Large spaces challenge us not only through the problem of diagonalizing large matrices, but also by requiring us to 
construct a good effective interaction. The bare nucleon-nucleon interaction needs to be modified in truncated spaces 
(this is an issue in the QRPA as well, though a less serious one). Currently, effective interactions are built through a 
combination of perturbation theory, phenomenology, and painstaking fitting. The last of these, in particular, becomes 
increasingly difficult when millions of matrix elements are required. 

Related to the problem of the effective interaction is the renormalization of transition operators. Though the 
problem of the effective Gamow- Teller operator which enters directly into 2v(3(5 decay, has drawn some attention, 
very little work has been done on the renormalization of the two-body operators that govern Qv(3(3 decay. Shell model 
calculations won't be truly reliable until they address this issue, which is connected with deficiencies in the wave 
function caused by neglect of single-particle levels far from the Fermi surface. Ref. [l^ suggests that significant 
improvement on the state of the art will be difficult but not impossible in the coming years. 



Constraining Matrix Elements with Other Observables 

The more observables a calculation can reproduce, the more trustworthy it becomes. And if the underlying model 
contains some free parameters, these observables can fix them. The renormalization of free parameters can make up 
for deficiencies in the model, reducing differences between, e.g., the QRPA and RQRPA once the parameters of both 
have been fit to relevant data. The more closely an observable resembles {)v(3(3 decay, the more relevant it is. 

Gamow- Teller distributions, both in the /3~ and (3^ directions, enter indirectly into both kinds of (3(3 decay, and 
are measurable through (p, n) and analogous nucleon exchange reactions. Ref.ll^l is particularly careful to reproduce 
those transitions as well as possible. Pion double charge exchange, in which a -k^ enters and a 7r~ leaves, involves 
the transformation of two neutrons into two protons, like (3(i decay, but the nuclear operators responsible aren't the 
same in the two cases. Perhaps the most relevant quantity for calibrating calculations of 0!//3/3 decay is 2v(3(5 decay, 
which has now been measured in 10 different nuclei. 

Two recent papers have tried to use 2u(3(3 decay to fix the strength of np pairing in QRPA-based calculations. Ref. 
[t^ used it only for the J'^ = l"*" channel relevant for 2v(3(3 decay, leaving the np pairing strength unrenormalized in 
other channels. By contrast, Ref Ijii renormalized the strength in all channels by the same amount. The results of 
the two procedures were different; the former reference found that the Qv(}(3 matrix elements depended significantly 
on the theoretical approach, while the later one found almost no dependence on model-space size, on the form of the 
nucleon-nucleon interaction, or on whether the QRPA or RQRPA was used. The authors argued that fixing the np 
pairing strength to 2v(3(3 rates essentially eliminates uncertainty associated with variations in QRPA calculations of 
01/(3(3 rates, though they left open the question of how close to reality the calculated rates were. 

The conclusion of [t^ is supported by the work "s^l where the 2v(3(3 was also used to fix the relevant parameter and 
essentially no difference between QRPA, RQRPA, and the new variant develo ped there was found. Moreover, when 
the fact that higher order terms are not included in [s^l, but are included in fa reduction by ~ 30%) are taken 
into account, these two calculations agree quite well. Another case is the work |81| which uses yet another variant of 
the theory and seemingly disagrees drastically with other calculations for the case of ^^Ge Qv(3(j decay. However, in 
an earlier work using the same method and parameters [s^, it is clear that the rate of the 2v(3(5 decay is incorrect 
by a large factor. When the relevant parameters are adjusted to give the correct 2v(3(3 rate, the result again agrees 
quite well with |79j . 
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While this question remains open, it is clear that only calculations of the OvP/S nuclear matrix elements that also 
describe other relevant nuclear properties should be included in the estimate of the uncertainty. 

Reducing the Uncertainty 

What can be done to improve the situation? In the near term, improvements can be made in both QRPA-based and 
shell-model calculations. First, existing calculations should be reexamined to check for consistency. One important 
issue is the proper value of the axial- vector coupling constant gA, which is often set to 1 (versus its real value of 
1.26) in calculations of /3 decay and 2^(3(3 decay to account for the observed quenching of low-energy Gamow- Teller 
strength. What value should one use for 0^(3(3 decay, which goes through intermediate states of all multipolarity, not 
just 1+? Some authors use = 1, some qa = 1.26, and some = 1 for the 1+ multipole and 1.26 for the others. 
(Often authors don't reveal their prescriptions.) The second of these choices appears inconsistent with the treatment 
of 2^(3(3 decay. Since the square of gA enters the matrix element, this issue is not totally trivial. The striking results of 
Ref. 79] suggest that an inconsistent treatment is responsible for some of the spread in the calculated nuclear matrix 
elements. More and better charge-exchange experiments would help teach us whether higher-multipole strength is 
also quenched. 

Next, the various versions of the QRPA should be tested against exact solutions in a solvable model that is as 
realistic as possible. The most realistic used so far are the S'0(5)-based model used to study the QRPA and RQRPA 
for Fermi 21/(3(3 decay a two- level version of that model used in 76] for the QRPA in Fermi 2i//3/3 and 01^/3(3 decay, 
and an S'0(8)-based model used to test the QRPA and RQRPA for both Fermi and Gamow- Teller 2i//3/3 decay in Ref. 
[s^ ]. It should be possible to extend the 5*0(8) model to several sets of levels and develop techniques for evaluating 
Qv[3(3 matrix elements in it. All these models, however, leave out spin-orbit splitting, which weakens the collectivity 
of np pairing. Using these models should help to understand the virtues and deficiencies of various QRPA extensions. 

Along the same lines, we will need to understand the extent to which such methods can reproduce other observables, 
and their sensitivity to remaining uncertainties in their parameters. A careful study of the first issue was made in 
These efforts must be extended. The work is painstaking, and perhaps not as much fun as concocting still 
more variations of the QRPA, but it is crucial if we are to reduce theoretical uncertainty. Self-consistent Skyrme 
HFB+QRPA, applied to single-/? decay in [^] and Gamow- Teller resonances in may be helpful here; it offers a 
more general framework, in principal anyway, for addressing the variability of calculated matrix elements. Solvable 
models can be useful here too, because they can sometimes supply synthetic data to which parameters can be adjusted 
(as in 76]). 

The best existing shell-model calculation produces smaller matrix elements than most QRPA calculations. Com- 
puter speed and memory is now at the point where the state of the shell-model art can be improved. The calculation 
of the (3(3 decay of ''^Ge in |j2J]used the /5/2P3/2P1/259/2 model space, allowing up to 8 particles (out of a maximum 
of 14) into the gg/2 level. Nowadays, with the help of the factorization method (see 0,123 )i S''^ accurate approx- 
imation to full shell-model calculations, we should be able to fully occupy the gg/2 level, and perhaps include the 
(77/2 and f'j/2 levels (though those complicate things by introducing spurious center-of-mass motion). In addition, one 
can try through diagrammatic perturbation theory to construct effective Qvl3(3 operators for the model space that 
are consistent with the effective interaction. Though perturbation theory has convergence problems, the procedure 
should at least give us an idea of the uncertainty in the final answers, perhaps also indicating whether result obtained 
from the "bare" operators is too large or too small. Research into effective operators has been revived in recent 
years jSd] and we can hope to improve on diagrammatic perturbation theory. One minor source of uncertainty con- 
nected with renormalization (which also affects the QRPA) is short-range two-nucleon correlations, currently treated 
phenomenologically, following [9(1 ]. 

In short, much can be done and we would be well served by coordinated attacks on these problems. There are 
relatively few theorists working in j3[3 decay, and their efforts have been fragmented. More collaborations, postdoctoral 
and Ph.D projects, meetings, etc., would make progress faster. There is reason to be hopeful that the uncertainty will 
be reduced. The shell-model matrix element may be too small because it does not include any particles outside the 
fp + f/g/2-shell. These particles, as shown by QRPA calculations, make the matrix element larger. We suspect that 
the results of a better shell- model calculation will be closer than the best current one to the QRPA results and that, 
as noted above, the spread in those results can be reduced. Finally, other nuclei may be more amenable to a good 
shell-model calculation than Ge. ^^^Xe has 82 neutrons (a magic number) making it a particularly good candidate. 
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IV. COSMOLOGY AND NEUTRINO MASS 

Neutrinos play an important role in astrophysics and cosmology (see also section HID). In cosmology, relic neutrinos 
may constitute an important fraction of the hot dark matter (HDM) influencing the evolution of large scale structures 
(LSS)^3|- The imprint of neutrino HDM on LSS evolution is quite distinct from other dark matter candidates such as 
supersymmetric particles, which act as non-relativistic or Cold Dark Matter (CDM). Cosmological models of structure 
formation strongly depend on the relative amounts of CDM and HDM in the universe, hence a determination of the 
neutrino contribution Qi, to the total dark matter content of the universe is important for our understanding of 
structure formation 92]. Such an investigation is a strong motivation for the next-generation terrestrial neutrino mass 
experiments. For example, a single /3-decay experiment with a neutrino mass sensitivity of 0.2 eV (90%CL.) will be 
sensitive to a neutrino contribution = 0.015. Such measurements would either significantly constrain or fix the role 
of HDM in structure formation by a model independent method. 

Two areas of astrophysics where an understanding of neutrino mass is important are supernova dynamics and 
cosmic ray physics. In the area of cosmic rays, a new model has been proposed which aims at explaining the origin 
of ultra high energy (UHE) cosmic rays. This so-called Z- burst model would require relic neutrino masses within the 
sensitivity range of next generation neutrino mass experimmcnts. 

Cosmological studies and neutrinos 

The contribution fli, of relic neutrinos to the total density Q can be probed by a combination of precise measurements 
of the temperature fluctuations of the cosmic microwave background radiation and with measurements of the matter 
fluctuation at large scales by high statistics galaxy redshift surveys. The results for the mass density of neutrinos can 
then be transformed into neutrino mass results based on the calculated relic neutrino density of 112j/'s/cm^ in the 
framework of the canonical standard cosmological model. (Detailed information on these neutrinos and cosmology 
can be found in the Astrophysics and Cosmology working group report.) 

Results from cosmological studies 

Early cosmological bounds on Emj^, based on pre-WMAP/SDSS/ 2dFGRS data, gave upper limits in the range of 
3-6 eV, comparable to the present laboratory limit from tritium /3-decay. Recent cosmological studies are based on 
different combinations of the high quality data from WMAP, small scale high resolution CMBR experiments, and the 
two large galaxy redhift surveys 2dFGRS and SDSS. In some cases additional structure information from Lyman-a 
data or X-ray luminosities from clusters have been added. 

Table IIIII presents a summary of either upper limits or best values for neutrino masses derived from recent cos- 
mological studies. The table shows the considerable spread in the results published recently. This is due to several 
generic difficulties associated with cosmological studies. First, these studies suffer from the problem of parameter de- 
generacy 97]. Different combinations of cosmological parameter can describe the LSS and CMBR data equally well, 
so additional information is required to break the degeneracy. The errors associated with these input parameters also 
imply that the v-mass results crucially depend on the priors for these parameters. Different priors for cosmological 
parameters result in limits for the sum of neutrino masses Em^ which differ by factors of 2 or more. 

The strong dependence of cosmological neutrino mass results can be illustrated by comparing the strongest upper 
limit on reported in the literature, the WMAP upper limit [s^ of m,y < 0.27 eV (95 % CL.), with tentative evidence 
for non-zero neutrino masses reported by Allan et al. Both analyses use substantially identical sets of input 

parameters, most notably the WMAP and 2dFGRS data. However, while the WMAP authors^^ use additional 
Lyman-a data for their analysis, Allan et al. use X-ray luminosity functions (XLF) of galaxy clusters obtained 
with the orbiting Chandra X-ray telescope. This small change of input data transforms an upper limit into evidence 
for non-zero masses with Em^=0.56eV as best fit value. It is interesting to note that both the use of Lyman-a 
data as well as the use of the measured XLF have been criticized in the literature. This clearly underlines that 
cosmological studies, as impressive as they are with regard to the determination of cosmological parameters, still yield 
model- dependent results for neutrino masses. 

Further concerns with cosmological neutrino mass studies are associated with systematic errors. LSS data from 
galaxy surveys suffer from the problem of biasing (i.e. to what extent does the galaxy distributing trace the distribution 
of cold dark matter and dark baryons), possible redshift-spacc distortions and selection effects. A model-independent 
input of the value of the neutrino mass to cosmological studies would thus be especially important for future high 
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author 


WMAP 




SDSS 


2dF 


other data 


Em, [eV] 


Bar'03 [93] 


X 


X 


X 


X 


h(HST) 


<0.75 


Teg'03 [94] 


X 


X 


X 




SNIa 


<1.7 


ASB'03 f95] 


X 


X 




X 


XLF 


= 0.36-1.03 


WMAP [32] 


X 


X 




X 


Lya , h (HST) 


<0.7 


Bla'03 [96] 


X 






X 


— 1 


= 2.4 


Han'03 [97] 


X 


X 




X 


h (HST), SNIa 


<1.01 


Han'03 [97] 


X 


X 




X 




< 1.2 


Han'03 [97] 


X 






X 




<2.12 



TABLE III: Survey of neutrino mass results obtained from the most recent cosmological studies. For each study the specific set 
of input data is listed: WMAP angular and/or polarisation data, high resolution CMBR experiments like CBI etc, LSS data 
from the SDSS and the 2dFGRS galaxy surveys, and other data. 

precision cosmological studies. A laboratory measurement of nii, could break the existing degeneracies between fi^ 
and other cosmological parameters, and thus help to provide a better picture of large scale structure evolution. 

Beyond the 'concordance' model 

The results of the cosmological studies presented above in table lllll have been obtained within the 'canonical standard 
cosmological model'. In the following we briefly list studies which are outside the present so-called 'concordance' flat 
ACDM models. 

In jo^ the authors argue that cosmological data can be fitted equally well in the framework of an Einstein-de Sitter 
universe with with zero cosmological constant A, albeit at the expense of requiring a very low value for the Hubble 
constant of Hq ~ 46 km/s/Mpc. The authors claim that CMB and LSS data seem to imply the existence of non-cold 
dark matter component. 

I n addition, several authors have speculated about the interesting coincidence between the smallness of neutrino 
masses (m^ < 1 eV) and the deduced vacuum density pv ~ ( 10"'^ eV )^ responsible for the cosmological dark energy A. 
In this framework, the masses of the active i^-species vary like the inverse of the I'-density. Accordingly, the neutrino 
mass could be of order of 1 eV today. 

Another 'coincidence' problem (J^a — ^m) has been used by M Tegmark et al. [9^ to investigate the possibility 
that mjy is a stochastic variable, which is randomized during inflation. 

Yet another possible modiflcation to the neutrino sector of the 'concordance' model has been brought up by the 
authors of [lOQ j. They assume that neutrinos have small extra scalar or pseudoscalar interactions. Should neutrinos 
possess even only tiny couplings of the order of 10^^ to hypothetical scalar bosons, the neutrino density il^ could be 
affected strongly by annihilation processes of the type v v *^ (j) (j). In this scenario neutrinos would not decouple from 
matter at T = lMeV, but would stay in thermal equilibrium until much later times (T = leV), which would inhibit 
free streaming. 

Relic neutrinos have not been detected at present, and will likely not be detected in the nearer future due to their 
very low energy in the /icV range. For this reason, the 'ncutrinolcss universe' scenario can only be tested by comparing 
the model-independent neutrino mass result from single and double beta decay experiments along with cosmological 
data. 

Future Perspectives 

At present, the cosmological studies of m^ are still model-dependent due to systematic effects such as biasing, 
parameter degeneracy, possible selection effects, possible contributions from non-linear effects and the strong influence 
of priors on the neutrino mass results. Future high precision studies aim at strongly reducing these systematic effects. 

In the field of CMBR experiments, the new promising technique of stud ying distortions of the CMBR temperature 
and polarization maps induced by gravitational lensing has been proposed lOl'. This technique is sensitive to changes 
of structure evolution at late times due to massive neutrinos and thus could break the degene racy between neutrino 
mass, equation-of-statc of the dark energy and the primordial power spectrum amplitude. In |lOll | the sensitivity of 
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this method is estimated to come down to Emi^=0.3 eV. This is of the same order as what is expected in |l02l| for the 
i/-mass sensitivity (m^=0.14 eV) from the Planck sateUite scheduled to start operations in 2008. 

These investigations will be complemented by future deep galaxy surveys extending out to redshift parameters z=2 
(DEEP2, VLT-Virmos) as well as dedicated studies of lensing effects on galaxy clusters (LSST, Large Synoptic Survey 
Telescope) . 

V. EXPERIMENTAL PROSPECTS FOR f3 DECAY 
Motivation for absolute measurements 

The SM of particle physics describes present experimental data up to the electroweak scale, but does not explain 
the observed pattern of the fermion masses or the mixing among the fermion generations. In addition it explicitly 
assumes neutrinos are massless and offers no explanation for the observed i^-masses and i/-mixing. 

There are many theories beyond the Standard Model, which explore the origins of neutrino masses and mixing. 
In these theories, which often work within the framework of Supersymmetry, neutrinos naturally acquire mass. A 
large group of models makes use of the so-called see-saw effect to generate neutrino masses. Other classes of theories 
are based on completely different possible origins of neutrino masses, such as radiative corrections arising from an 
extended Higgs sector. As neutrino masses are much smaller than the masses of the other fermions, the knowledge 
of the absolute values of neutrino masses is crucial for our understanding of the fermion masses in general. Recently 
it has been pointed out that the absolute mass scale of neutrinos may be even more significant and straightforward 
for the fundamental theory of fermion masses than the determination of the neutrino mixing angles and CP-violating 
phases j91j. It will likely be the absolute mass scale of neutrinos which will determine the scale of new physics. 

Introduction (related to kinematical experimental tectiniques) 

Two kinematical experimental techniques that probe neutrino mass directly, are precise observations of decay 
kinematics in nuclear or particle decay processes and the utilization of time of flight measurements in the detection of 
supernova neutrinos incident on terrestrial neutrino detectors. The decay measurements, which are based on purely 
kinematical observables, have essentially no reliance on theoretical assumptions about neutrino properties and hence 
offer a completely model independent probe of absolute neutrino mass. 

Examples of decay measurements include studies of f3 -decay spectra shapes, measurements of muon momentum in 
pion decay, and invariant mass studies of multi-particle semileptonic decays of the r. However, the recent measurement 
of large mixing angles in the lepton sector gives a clear advantage to (3 -decay spectra shape measurements, since these 
experiments probe all three neutrino mass eigenstates and are approaching the sub-eV range os sensitivity. These 
next generation (3 -decay measurements will complement other laboratory and cosmological methods to investigate 
neutrino masses. The combination and comparison of results from /3-decay, 0v[3(3 and cosmological studies will be 
essential for our understanding of the role of neutrinos in our physical world, both at the Micro- and Macro world. 

In contrast to 0i>(3P experiments, kinematic investigations of the neutrino mass do not rely on further assumptions 
on the neutrino mass type, Majorana or Dirac. Such kinematic experiments can be classified into two categories 
both making use of the relativistic energy momentum relation = p^c^ + m^c* as well as of energy and momentum 
conservation. 

Neutrino time-of- flight studies 

The narrow time signal of a supernova (SN) neutrino burst of less than 10 s in combination with the very long- 
baseline between source and detector of several kpc would allow the investigation of small ToF effects resulting from 
small z/-masses. Supernova ToF studies are based on the observation of the energy-dependent time delay of massive 
neutrinos relative to massless neutrinos. This method provides an experimental sensitivity for the rest masses of i/g, 
I'^and VrOi a few tens of eV. This sensitivity can be pushed into the few-eV range, if additional assumptions on the 
time evolution of the z^-burst are being made. In this case the v-va.ass sensitivity becomes model-dependent, however. 

Recently, new methods have been proposed which do not rely on details of the i^-burst timing which might yield 
sensitivity in the few-eV range for and Vr- The two most promising techniques arc: a) the measurement of the 
abrupt termination of the SN-i' signal due to the early formation of a black hole (161- .162j and b) the correlation 
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of SN i/-signals with independent signals from gravitational wave experiments |l63l | . Except for exceedingly close by 
supernova, < 1 kpc , the sensitivity from the first technique to neutrino mass would still only be at the level of a few 
eV. Estimates for the sensitivity to the second method are also expected to also be only at the eV level. 

If future measurements were to reveal a SN i^-pulse which is terminated abruptly after a few seconds, the infor- 
mation on the i/-mass provided by next-generation terrestrial experiments would likely be used to help further our 
understanding of supernovae dynamics. 



Weak particle decay processes 

The investigation of the kinematics of weak decays is based on the measurement of the charged decay products 
of weak decays. For the masses of and Vr the measurement of pion decays into muons and at PSI and the 
investigation of r-decays into 5 pions and at LEP have yielded the upper limits: 

TO(t/^) < 190 keV at 90 % confidence 
m{v-r) < 18.2 MeV at 95 % confidence 

Both limits are much larger than those attainable in beta decay spectra studies. 

P decay 

The technique for detecting neutrino mass in beta decay is essentially to search for a distortion in the shape of the 
beta spectrum in the endpoint energy region. The most sensitive searches for the electron neutrino mass up to now 
are based on the investigation of the electron spectrum of tritium /? decay 

Hc+ + e- +De . (34) 
The electron energy spectrum of tritium (3 decay for a neutrino with mass is given by 

^ = C X F{Z, E)pEiEo - E)[{Eo - Ef - ml]ie{Eo -E-m,), (35) 

where E denotes the electron energy, p is the electron momentum, E^ corresponds to the total decay energy, F{Z^ E) 
is the Fermi function, taking into account the Coulomb interaction of the outgoing electron in the final state, the 
stepfunction 0(i?o — E — rriu) ensures energy conservation, and C is given by 

C^Gl-ArCosHcWl^" . (36) 

Here Gf is the Fermi constant, 9c is the Cabibbo angle, rue the mass of the electron and M is the nuclear matrix 
element. As both M and F{Z, E) are independent of m^, the dependence of the spectral shape on rrii, is given by the 
phase space factor only. In addition, the bound on the neutrino mass from tritium (3 decay is independent of whether 
the electron neutrino is a Majorana or a Dirac particle. 

The signature of an electron neutrino with a mass of m(i^e )=leV is shown in Fig. |31 in comparison with the 
undistorted (3 spectrum of a massless v,, . The spectral distortion is statistically significant only in a region close to 
the /3 endpoint. This is due to the rapidly rising count rate below the endpoint dN/dE ex (i?o ~ E)"^ . Therefore, only 
a very narrow region close to the endpoint Eq is analyzed. Note that in fitting the measured shape to a calculated 
shape the functional form depends on and is not defined for < 0. 

One immediate difficulty apparent in the figure is that there are very few decays in the region of interest, only 
2 • 10^^^ of the total rate is in the last eV. A related problem is that one must carefully eliminate or minimize 
backgrounds. In addition subtle effects such as instrumental resolution, energy loss of the electrons, and atomic or 
molecular physics excitations during the decay can cause shape distortions that are of similar size to the effect of 
nonzero neutrino mass but of opposite sign. Hence tritium (3 decay experiments with high neutrino mass sensitivity 
require a huge luminosity combined with very high energy resolution. Furthermore, a precise determination of a value 
or limit for requires complete and accurate understanding of all systematic effects that can alter the shape of the 
spectrum. 

Apart from offering a low endpoint energy Eq and a moderate half life of 12.3 y, tritium has further advantages as 
(3 emitter in v mass investigations: 
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1. the hydrogen isotope tritium and its daughter, the '^He^ ion, have a simple electron shell configuration. Atomic 
corrections for the (3 decaying atom -or molecule- and corrections due to the interaction of the outgoing /3-electron 
with the tritium source can be calculated in a simple and straightforward manner 

2. The tritium [3 decay is a super-allowed nuclear transition. Therefore, no corrections from the nuclear transition 
matrix elements M have to be taken into account. 

The combination of all these features makes tritium an almost ideal (3 emitter for neutrino mass investigations. 



Current tritium /3-decay results 



The Mainz and Troitsk groups have set the most precise limits on the electron antineutrino mass. Both experiments 
utilize novel magnetic solenoidal retarding electrostatic spectrometers which measure an integral beta spectrum, 
integrating all energies above the acceptance energy of the spectrometer. In their measurements, the Mainz group 
utilized a frozen molecular tritium source. Their result Il65l is: 



m2_^ = -1.2±2.2±2.1 (37) 

which yields a limit of: 

m^^ < 2.2 eV {95%CL). (38) 



This result is based on data that has passed several systematic and consistency checks. The Troitsk grouD |16fil ll 
developed a gaseous molecular tritium source and has also published a limit similar to that of the Mainz group of 

m2_^ = -2.3±2.5±2.0 ey^ (39) 

with a limit of: 

m^^ < 2.1 eV (95%CL). (40) 
However, they must include a not well understood step function near the endpoint in order to produce such a limit. 



Next generation experiments 



The K.Arlsruhe THItium Neutrino project (KATRIN) experiment 



The KArlsruhe TRItium Neutrino project (KATRIN) experiment is a ne xt-generation tritium /3-decay experiment 
designed to measure the mass of the neutrino with sub-eV sensitivitv |l68j . KATRIN utilizes a windowless gaseous 
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molecular tritium source and a 10 m diameter magnetic solenoidial retarding electrostatic type spectrometer. The 
experiment will be constructed adjacent to the Tritium Laboratory Karlsruhue at the Forschungszentrum Karlsruhe 
(FZK). The collaboration membership includes a significant number of collaborators from the most recent successful 
tritium /3-decay experiments, including Mainz and Troitsk, as well as the Washington group which developed the 
original LANL gaseous based tritium source. The current schedule calls for for KATRIN to become operational 
during 2008. The experiment expects to achieve an order of magnitude more sensitivity than the best /3-decay 
measurements carried out to date. (Note, since the shape effects are proportional to this implies a factor of 
100 increase in sensitivity.) Based on the current design, KATRIN expects after three years of running to reach a 
sensitivity to neutrino mass of 0.20 eV (90% CL) and would hence be able to observe a neutrino mass with a mass of 
0.35eV at the 5 sigma significance level. The experiment has recently received a firm funding commitment from the 
German funding ministry. 



NEXTEX (Neutrino Experiment at Texas) 



The Neutrino Experiment at TEXas is an experiment that has been under construction for more than 10 years 
with the goal of measuring the electron antineutrino mass to 0.5 eV (lower limit at 3(t). Collaborating institutions 
are The University of Texas at Austin, Nebraska, Brandeis, Michigan Technological Institute, Pomona College, and 
Southwestern University. The experiment uses the end point method for gaseous tritium decay, and the apparatus con- 
sists of three serial electrostatic differential analyzers. Presently the experiment is funded by the NSF to demonstrate 
proof of principle. Construction completion and final approval by the site committee is expected by the collaborators 
within 1 year. The experiment will run for 3 years and the anticipated funding request would be $3.9M. 



Other approaches to 13 decay 

A different approach to directly measure the electron neutrino mass is the use of cryogenic bolometers. In this case, 
the /? source can be identical to the /3-electron spectrometer. Thi s new technique has been applied to the isotope 
Re, which has a 7 times lower endpoint energy than tritium.l69ill7nl|. Current microcalorim eters reac h an energy 
resolution of /S.E ~ 5 eV for short-term measurements and yield an upper limit of m(i/e) < 15 eV [TfiilrTfij . To further 
improve the statistical accuracy, the principle of integration of active source and detector requires the operation of 
large arrays of microcalorimeters. The sensitivity on the neutrino mass in the next 5-7 years is expected to reach 
below the few eV levelll69t. 



VI. EXPERIMENTAL PROSPECTS FOR (5(3 



If an experiment observes Qv(3(3 it will ha ve pr ofound physics implications. Such an extraordinary claim will require 
extraordinary evidence. The recent claim |l03j for an observation of Ov(3(3 has been controversial (See discussion 
below). Also previous "false peaks" in (3(3 spectra have appeared near a Ov(3(3 endpoint energy (see discussion in 
|l04) . page 273). One must ask the question: What evidence is required to convincingly demonstrate that Ovf3(3 has 
been observed? Low-statistical-significance peaks (w 2a) have faded with additional data, so one must require strong 
statistical significance (perhaps 5cr). (See Fig. 2|) This will require a large signal-to-noise ratio that will most likely 
be accomplished by an ultra-low-background experiment whose source is its detector. Such experiments are usually 
calorimetric and provide little information beyond just the energy measurement. 

How does an experiment demonstrate that an observed peak is actually due to (3(3 decay and not some unknown 
radioactivity? Additional information beyond just an energy measurement may be required. For example, although 
there is some uncertainty associated with the matrix elements, it is not so large that a comparison of measured 
rates in two different isotopes could not be used to demonstrate consistency with the Majorana-neutrino hypothesis. 
Alternatively, experiments that provide an additional handle on the signal, for example by measuring a variety of 
kinematical variables, demonstrating that 2 electrons are present in the final state, observing the 7 rays associated 
with an excited state, or identifying the daughter nucleus, may lend further credibility to a claim. Experiments 
that provide this extra handle may require a significantly more complicated apparatus and therefore face additional 
challenges. 

The exciting aspect of (3(3 research today is that many proposed experiments intend to reach a Majorana mass 
sensitivity of y^Sm^^. Several different isotopes and experimental techniques are being pursued actively and many 
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FIG. 4: This contour plot shows the half life, in units of 10^^ y, for a peak of 5a significance for a given exposure and background. 
The KDHK point is shown. 

TABLE IV: A summary of the recent Oi'/3/3 results. The {m/sfs) limits are those deduced by the authors. All limits are at 90% 
confidence level unless otherwise indicated. The columns providing the exposure and background are based on arithmetic done 
by the authors of this paper, who take responsibility for any errors in interpreting data from the original sources. 



Isotope Exposure 


Background Half-Life 






(kmole-y) 


(counts) 


Limit (y) 


(meV) 


«Ca 


5 X 10"^ 





> 1.4 X 10^2 


< 7200 - 44700fl05l 


^"^Ge 


0.467 


21 


> 1.9 X lO^'' 


< 350|12^ 


^•^Ge 


0.117 


3.5 


> 1.6 X 10^^ 


< 330 - 1350fl07l 


^«Ge 


0.943 


61 


= 1.2 X 10^^ 


= 440[103] 


^'Se 


7 X 10"^ 





> 2.7 X 10^2(68%) 


< 5000fl08l 


io»Mo 


5 X 10~* 


4 


> 5.5 X 10^2 


< 2100fl09l 


"«Cd 


1 X 10~^ 


14 


> 1.7 X 10^^ 


< 1700fll0l 


Te 


Geochem. 


NA 


> 7.7 X lO^'' 


< 1100 - 1500flll1 


130 


0.025 


5 


> 5.5 X 10^^ 


< 370 - 1900fll2l 


"6Xe 


7 X 10"^ 


16 


> 4.4 X 10^^ 


< 1800 - 5200[1131 


isoNd 


6 X 10"^ 





> 1.2 X 10^1 


< 3000fll4l 



of the programs look viable. In this section we describe the current situation in experimental P(3 decay . 



Results to date 

Table IIVI lists the recent Oi^/3/3 results. The best limits to date come from the enriched Ge experiments. The 
two experiments had comparable results although the Heidelberg-Moscow result was marginally better. The 
limits near 2 x 10^^ y results in a {mpfj) li mit n ear 300 meV, with an uncertainty of about a factor of 3 because of 
the uncertainty in |Moi/|. One recent paper [lisj performed a joint analysis of the two experiments and found T^j2 

> 2.5 X 10^5 y. 

Most of the results listed in Table Hvl are at least a few years old. The obvious exceptions to this are the Te and 
Cd results. CUORICINO continues to collect data. 
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FIG. 5: The spectrum from the Heidelberg-Moscow experiment upon which the claim for 0/^/3/3 is based. The data in the two 
panels are identical. The lower panel has a Gaussian curve to indicate the strength of the claimed QuPP peak. 



A claim for the observation of Oz//3/3 



In early 2002, a claim for the observation of Qv(3l3 was p ublished ( Klap dor-Kleingrothaus et aL2002a). The paper 
made a p oor c ase for the claim and drew strong criticism 1 11 fj . 111?! Ill8l |. The initia l response to the criticism was 
emotional |ll9j . In addition, one of the original co-authors wrote a separate replv |l2(il | that mostly defended the claim 
yet acknowledged some significant difficulty with the analysis. This author's name doesn't appear on later papers. 
More rece ntly, however, supporting evidence for the claim has been presented and we recommend the reader study 
Ref. jl2ll | for a good discussion of the initial evidence and Ref. |l03| for the most recent data analysis. Importantly, 
this later paper includes additional data and therefore an increase in the statistics of the claim. In this subsection we 
summarize the current situation. (We use the shorthand KDHK to refer to the collection of papers supporting the 
claim.) 

Figure 13 shows t he s pectrum corresponding to 71.7 kg-y of data from the Heidelberg-Moscow experiment between 
2000 and 2060 keV |l03| . This spectrum is shown here to assist the casual reader in understanding the issues. However, 
the critical reader is encouraged to read the papers listed in the references as the authors analyze several variations 
of this data using different techniques. The fit about the expected Oj//3/3 peak energy yields 28.75 ± 6.86 counts 
assigned to Qu(3[3. The paper claims a significance of approxmately 4(t for the peak, where the precise significance 
value depends on the details of the analysis. The corresponding best-fit lifetime, Tjy'j = 1-19 x 10^^ vears^l03j. leads 
to a (m/3/3) of 440 meV with the matrix element calculation of Ref. |l22l | chosen by the authors. 

In the region between 2000 and 2100 keV, the KDHK analysis of 2002 found a total of 7 peaks. Four of these 
were attributed to '^^^B\ (2011, 2017, 2022, 2053 keV), one was attributed to OvjSP decay (2039 keV), and two were 
unidentified (2066 and 2075 keV). The KDHK analysis of 2004 does not d iscus s the spectrum above 2060 in detail. 
An additional possible feature may also be present near 2030 keV. A studv [l46j | comparing simulation to calibration 
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with ^-'^^Bi demonstrates that if the location of the Bi is known, the spectrum can be calculated. Furthermore, the 
relative strengths of the strong Bi lines at 609, 1764 and 2204 keV can be used to determine the location of the activity. 
Because the results of summing effects depend on the proximity of the activity, its location is critical for the simulation 
of the weak peaks near the OfP/S endpoint. The study also shows that the spectrum can't be naively estimated, as 
was done in Ref. |ll7l |. In fact, Table VII in Klapdor-Kleingrothaus (2002c) finds, even with a careful simulation, 
that the expected strengths of the ^^^Bi peaks in the 2000-2100 keV region are not predicted well by scaling to the 
strong peaks. That is, the measured intensities of the weak peaks are difficult to simulate without knowing the exact 
location of the activity. Furthermore, the deduced strengths of the weak lines are more intense than expected by 
scaling from the strong peaks, even though the activity location is chosen to best describe the relative intensities of 
the strong peaks. 

Double-beta decay produces two electrons that have a short range in solid Ge. Therefore, the energy deposit is 
inherently localized. Background process, such as the 7 rays from ^^^Bi, tend to produce multiple energy deposits. The 
pulse waveform can be analyzed to distinguish single site events (SSE) from multiple site events. Such an analysis by 
KDHK [iO^ I12.1I | tends to indicate that the Bi lines and the unidentified lines behave as multiple site events, whereas 
the QiyfUP candidate events behave as SSE. Note, however, that the statistics are still poor for the experimental lines 
and this conclusion has a large uncertainty. Nonetheless, this feature of the data is very intriguing and clearly a 
strength of the KDHK analysis. 

An analysis by Zdesenko et aL |ll5| points out the strong dependence of the result on the choice of the window width 
of the earlier 2002 analysis. The KDHK analysis argues that a small window is required because of the neighboring 
background lines. Even so, their Monte Carlo analysis shows that the result becomes less stable for small windows 
(see Fig. 9 in Ref. [mj). Zdesenko et al.^n^ also remind us that the significance of a signal is overesti mated when 
the regions used to estimate the background are comparable to the region used to determine the signal [l2^ . The 
report of Ref. (103j fits a wide region containing several peaks simultaneously after using a Bayesian procedure to 
identify the location of the peaks. 

The claim for OvjSf] decay was made by a fraction of the Heidelberg-Moscow collaboration. A separate group of 
the original collab oration presented their analysis of the data at the IV International Conference on Non- Accelerator 
New Phvsics [l25j . They indicate that the data can be separated into two distinct sets with different experimental 
conditions. One set includes events that are described as "underthreshold pulses" and one set that does not. Analysis 
of the two sets produce very different conclusions about the presence of the claimed peak. They conclude that the 
evidence is an experimental artifact and not a result of 0;//3/3decay. KDHK resDonds |lO,lj that these corrupt data 
were not included in their analysis. 

Traditionally, P(3 experiments have ignored systematic uncertainties in their analysis. Only recently with the start- 
up of high-statistics 2i//3/3 results has this situation begun to change. Historically, OvP/S results have always been 
quoted as upper limits based on low count rates. As a result, systematic uncertainties tended to be negligible in the 
final quoted values. With a claim of a positive result, however, the stakes are dramatically raised. It is clear that it 
is difficult to produce a convincing result when the signal counts are comparable to expected statistical fluctuations 
in the background. The further presence of nearby unidentified peaks makes the case even harder to prove. Although 
KDHK does discuss some systematic uncertainties qualitatively and indicates they are small (in the position of the 
OvfiP peak, and the expected peak width, for example), there is no consideration of an uncertainty associated with 
the background model. 

The next round of proposed 0z//3/3 experiments are designed to reach (5m^t,n and therefore will quickly confirm 
or repudiate this claim. This is fortunate since the feature near 2039 keV in the KDHK claim will likely require an 
experimental test. These experiments should provide a detailed listing of all identified systematic uncertainties and a 
quantified estimate of their size. Furthermore, because the stakes are very high and there will be many people who 
are biased, either for or against the KDHK claim, blind analyses should also become part of the experimental design. 



The Search for Decays to Excited States 

Searches for P(3 to excited states in the daughter atom have been p erform ed in a nu mber of isotopes but only 
observed in ^"^Mo (The experimental situation is reviewed by Barabash |l26j .') and ^^"Nd [l2^ . These experiments 
typically search for the 7 rays that characterize the excited states and therefore are not mode-specific searches. The 
interpretation therefore is that the measured rate (or limit) is for the 2u(3[3 mode. These data may be very useful to 
QRPA nuclear theory because the behavior of the nuclear matrix elemen ts with respect to gpp for the excited state 
decays is different than for transitions to the ground state |l28l Il29l Il30| . Thus, the excited state transitions probe 
different aspects of the theory and may provide insight into the physics of the matrix elements. 
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A further reason for interest i n decays t o the excited state, as mentioned earher, is the potential abihty to discover 
the process mediating t he decavllSll Il3^ . However, the decay rate to an excited state is 10-100 times smaher than 
rate to the ground state |l33l Furthermore the structure of the excited state in the daughter nucleus is not as 

well understood as the ground state, and this increases the relative uncertainty in the nuclear matrix element. 

The Search for Modes of Decay 

The f3~^ modes of decay have not received the attention of the modes because of the greatly reduced phase 
space and corresponding long half-lives. However, the rate is proportional to (m^/j)^ and its detection would provide 
additional matrix-element data. Furthermore, if the zero-neutrino mode were detected, it mig ht p rovide a handle on 
whether the decay is predominantly mediated by a light neutrino or by right-handed currents [l35j- 

Radiative neutrinoless double electron capture is a possible alternative to traditional neutrinoless double beta 
decaV [136j| . In this process, two electrons are captured from the atomic electron cloud and a radiated photon carries 
the full Q-value for the decay. A resonance condition can enhance the rate when the energy release is close to the 
2P-1S energy difference. In this case, high-Z, low-energy-release isotopes are favored (e.g. ^^^Sn). Unfortunately the 
mass differences for the candidate isotopes are not known precisely enough to accurately predict the overlap between 
the two energies. If a favorable overlap does exist, however, the sensitivity to {mpjs) might rival that of Oj^/3/?decay. 

Towards a 100-kg experiment 

The KDHK spectrum shows a feature very close to the Qvfip endpoint. This intriguing result will need to be 
confirmed or refuted experimentally. One can see the required operation parameters for a confirmation experiment 
from the KDHK result. One needs about 75 kg-y of exposure, and a background lower than about 0.5 counts/(kg y). 
Note that most of the proposals described above will all accomplish this very early on in their program if they meet 
their design goals. If instead one designs an experiment only to test the claim (not to provide a precise measurement 
of the T^j2) then a 100-kg experiment could provide the answer after a modest run time. 

If the KDHK result holds up, it will be a very exciting time for neutrino-mass research. A {mpp) near 400 meV 
means that /?-decay experiments and cosmology will be sensitive to the mass. As a result, one can certainly imagine 
a not-too-distant future in which we know the neutrino mass and its Majorana-Dirac character. Towards this goal, a 
precision measurement of (mfs^) will be required. To accomplish this, we will need more than one experiment, each 
with a half-life measurement accurate to 10-20%. At this level the uncertainty will be dominated by the matrix element 
uncertainty even if future calculations can be trusted to 50%. With two experiments utilizing different isotopes, one 
might disentangle the uncertainty in |Mq^|. 

Towards a 100-ton experiment 



The next generation of experiments hopes to be sensitive to v^^^^ltm- they fail to see 0i//3/3 at that level, the 
target for the succeeding generation of efforts will be This scale is an order of magnitude lower and hence 

will require two orders of magnitude more isotopic mass, approximately 100 tons of isotope. 

Proposed Experiments for /3/3 Decay 

The recent reviews by Elliott and Vogel 1] and Elliott and Engelj^i] describe the basics of experimental Oi^/3/3 decay 
in some detail. Therefore, we refer the reader to those articles and only summarize the status of the various projects. 
Table fvl lists the proposals. 

CANDLES 

The CANDLES collaboration has recently published the best hmit on OvP/S decay of 1.4 x 10^^ y in "^^Ca llOSl . Using 
the ELEGANTS VI detector, this experiment consisted of 6.66 kg of CaF2(Eu) crystals surrounded by Csl crystals, 
a layer of Cd, a layer of Pb, a layer of Cu, and a layer of LiH-loaded paraffin, all enclosed within an air-tight box. 
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TABLE V: A summary of the Oi//3/3 proposals. Background estimates were not available for all projects. The quantity of 
isotope includes the estimated efficiency for OufS/S. 
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liquid Xe 



This box was then surrounded by boron-loaded water tanks and situated underground at the Oto Cosmo Observatory. 
This measurement successfully demonstrated the use of these crystals for stti dies. 

An improved version of this crystal technology, the CANDLES-III detector Ts^, is presently being constructed with 
200 kg of CaF2 crystals. These crystals have better light transmission than the CaF2(Eu) crystals. This design uses 
sixty 10-cm'^ CaF2 crystals, which are immersed in liquid scintillator. The collaboration has also proposed a 3.2-t 
experiment that hopes to reach 100 meV for {mpp). 



COBRA 



The COBRA experiment [l40| uses CdZnTe or CdTe semiconductor crystals. These crystals have many of the 
advantages of Ge detectors but, in addition, operate at room temperature. Because the crystals contain Cd and 
Te, there are 7 00 and 0~^0^ isotopes contained. The final proposed configuration is for 64000 1-cm'^ crystals for 
a total mass of 370 kg. The collaboration has already obtained 30-keV resolution at 2.6 MeV with these detectors 
and has published initial /3/3-decay studies [l5^ . Background studies are the current focus of the efforts. Although 
it is tempting to focus on the naturally isotopic abundant ^^°Te for 01/00 decay, the presence of the higher Q-value 
^^^Cd creates a serious background from its 21/00 decay. Detectors enriched in ^^^Cd are probably required to reach 
45 meV. 



CUORE 



CUORE and CUORICINO are based on the technique of cryogenic detectors. When operated at low temperature, 
the absorbers of these detectors have a heat capacity so low that even the small energy released by a single radioactive 
decay event can be observed and measured by means of a suitable thermal sensor. With crystals of mass near to a 
kilogram, with NTD Ge (Neutron Transmutation Doped germanium) thermistors, an energy resolution similar to that 
of germanium diodes has been achieved. In addition, thermal detectors allow a wide choice of nuclei to be used for 
double beta decay searches. The experiment CUORICINO is located in the Gran Sasso underground laboratory and 
it is a prototype for CUORE (Cryogenic Underground Observatory for Rare Events). CUORICINO is an array of 44 
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crystals of Te02 each 5x5x5 cm and 18 crystals each 3x3x6 cm. With its mass of approximately 40 kg, CUORICINO 
is by far the most massive cryogenic set-up in operation. Due the large isotopic abundance (34%) of the double beta 
decay candidate ^^°Te, no isotopic enrichment is required, but two of the 3x3x6 cm crystals are enriched in ^-^^Te 
and two other in ^^*Te to investigate 2^/3/3. In only three months of operati on, C UORICINO has obtained a 90% c.l. 
limit on the lifetime against neutrinoless double beta decay of 7.5 x 10^'^ vr |ll2| . corresponding to an upper limit on 
the average neutrino mass ranging from 0.3 to 1.7 eV. This result rivals the best limits obtained from many years of 
searches for the double beta decay of ^^Ge. CUORE will consist of an array consisting of 25 columns of 10 planes of 
4 Te02 crystals each 5x5x5 cm of for a total of 1000 crystals with a mass around 760 kg. Each tower will therefore be 
similar to the single tower of CUORICINO, which consists of 13 planes. As far as time is concerned, CUORICINO 
is now running and its larger brother CUORE will be available in four years from the start of construction (likely 
summer of 2004). 

The present background of CUORICINO in the region of neutrinoless double beta decay (0.20+0.02 
counts/keV/kg/year) is in excellent agreement with the previously predicted value (0.22 counts/keV/kg/year). In the 
present measurement this background is mainly due to the surface contamination of copper and crystals and recently 
it has become understood how to reduce it by an order of magnitude by surface treatment. Taking into account that 
the structure of CUORE allows a large suppression of background by applying the anticoincidence method we can 
guarantee a conservative value of background of 0.01 counts/keV/kg/year. We believe that in the next four years we 
can achieve a further improvement in the energy resolution, in the radioactive contamination, and in the neutron and 
cosmic ray background. Thus, it is reasonable to predict a background of 0.001 counts/keV/kg/year for CUORE. 
As a conseque nce we believe that the CUORE sensitivity can be in the few tens of milliclcctronvolts for the average 
neutrino massll41l. 



DCBA 

As demonstrated by the history of discoveries of extremely rare events, magnetic tracking detectors have played 
important roles. Thus this collaboration bases its search for 0i>(3P e vents on a magnetic tracking detector. A momen- 
tum analyzer called the Drift Chamber Beta-ray Analyzer (DCBA) jl42| is a tracking detector operated in a uniform 
magnetic field of around 1 kC Various isotopes can be installed within the detector, if the source can be fabricated 
into a thin plate. Presently, the collaboration is considering ^^Se, ^""Mo, and ^^"Nd because of their high Q-values. 
A tracking region on each side of a source plate includes anode-, potential- and cathode- wires. The drift region is 
filled with 1-atm helium gas mixed with small amounts of a quench gas. A /3 ray emitted from a source plate makes 
a helical track in the region between anode and cathode wire planes. Anode signals are read out with Flash Analog 
to Digital Converter (FADC). The three-dimensional reconstruction of a helical track is available using data from 
the electron drift time (corresponding to the X-coordinate) , an anode wire position (Y) and the ratio of signals from 
both sides of an anode wire (Z). Momentum of each (3 ray is derived from the curvature of the track. Since financial 
support for DCBA has not been approved yet, the construction schedule is unknown. Electron tracks of 1 MeV were 
studied using internal conversion electrons from ^°^Bi, which was installed in a prototype called DCBA-T. We are 
making efforts to improve the Z-position resolution so as to obtain better energy resolution. Another developing item 
is to accommodate source plates, as much as possible, in a limited chamber volume. Research into cleaning the source 
material is also proceeding. 

The future DCBA experiment will consist of 40 modules. One module comprises a drift chamber of about 1.8 
m^ volume containing the source plates surrounded by a solenoid magnet with maximum field of 1.6 kC. In the 
chamber, 30 tracking regions cover 29 source plates. The total source-plate area is 25 m^ in each module. For a 
source-plate thickness of 60 mg/cm^, the source weight is 15 kg for each module. Therefore total source weight is 
600 kg, corresponding to about 6600 mol and 5400 mol for *^Se and ^'^°Mo, respectively, which are enriched to 90%. 
For a natural Nd source, about 200 mol of ^^°Nd will be installed. Assuming an efficiency of event detection of 0.3, a 
background rate of 1 event/module/year, and a measuring time of 5 years; the half- life sensitivities are approximately 
calculated to be 3 x lO^^ yr for ^^Se, 2 x lO^^ yr for ^°°Mo and 9 x lO^" yr for ^^ONd in natural Nd. If 90% enriched 
^^°Nd is available in the future by the method of Atomic Vapor Laser Isotope Separation, it is possible to obtain 
1 x 10^6 yr for isoNd. 
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EXO 

The Enriched Xenon Observatorv |l43j is being set-up to study the double beta decay of ^'^^Xe. EXO is a coUabo- 
rative effort currently involving groups from Caltech, Carleton University (Canada), Colorado State University, ITEP 
Moscow (Russia), Laurentian University (Canada), SLAC, Stanford University, University of Alabama, University of 
California Irvine, Universite de Montreal (Canada), Universite de Neuchatel and Universite de Yverdon (Switzerland). 
The collaboration aims to use up to 10 tons of isotopically enriched ^^^Xe to build a redundant and background-free 
detector using good energy resolution, pattern recognition and the identification of th e atomic species produced by 
the double-beta decay. A concise description of the project has been published in |l43l |. 

The final state atom tagging is possible because of the simple and well known atomic spectroscopy of Ba+ ions. Such 
spectroscopy has enabled the observation of individual ions illuminated with appropriate wavelengths since about 20 
years. The specific wavelengths needed to produce atomic fluorescence ensure extreme selectivity of this technique. 
Ba happens to be the atomic species produced in the double beta decay of ^^^Xe. In EXO the Xenon will be used as 
an active target in a Time Projection Chamber (TPC) either in liquid (LXe) or gas (GXe) phase. In the GXe case 
the laser beams would be steered to the location where a candidate decay has occurred. In the LXe case the Ba-ion 
candidate would be extracted and brought into an ion trap where the fluorescence would be observed. The possibility 
of observing the fluorescence of the Ba directly in the liquid is also being investigated by one of the EXO groups. 
While R&D is proceeding at different institutions for both liquid and gas phase TPC, a LXe TPC for 200 kg of ^^^Xe 
is being built as a prototype and as a first step towards the very large detector. 

EXO plans to test Majorana neutrino masses as small as 10 to 40 meV using this scheme |l43l | that should ensure 
extremely high background rejection power. This sensitivity covers neutrino masses derived from the atmospheric 
mass spitting. A degenerate or inversely hierarchical neutrino mass pattern can thus be tested with EXO. 

^•^^Xe is a particularly convenient isotope for a very large double beta decay experiment. It combines a large Q- 
value with ease (and low cost) enrichment, the absence of long lived radioactive isotopes, ease of purification and the 
possibility of transfer from one detector to another in the case new technology would become available. In addition 
purification can be achieved on-line so that more refined purification system may be introduced as they become 
needed. Finally, together with the enabling possibility of final state identification that is the hallmark of EXO, ^"^^Xe 
has the longest 2j//3/3half life among all high Q-value /3/3-unstable nuchdes (at least factor 6.5 longer than e.g. ''^Ge). 

With its 200 kg of source strength the EXO prototype detector will already represent the largest existing double 
beta decay experiment. Recently evidence for a Majorana neutrino mass of 0.441q'2o eV (error s giv en at Scr c.l.) 
has been claimed by a part of the Heidelberg component of the Heidelberg-Moscow experiment |l03l |. Should this 
claim turn out to be correct the EXO pr ototy pe expects to observe 43I3Q 0^^/3/3 decays per year when using the same 
matrix element calculation as reference |l03| . Assuming that funding continues to be available in a timely fashion, 
the prototype detector is expected to be commissioned at Stanford in the Winter 2004-5 and be transferred to WIPP 
in the Summer 2005. The development of the barium extraction and tagging will continue in parallel, together with 
the work on a GXe TPC. After initial operation of the prototype the technology for the large detector will be chosen 
and its design finalized. 



MOON 

MOON (Molybdenum Observatory Of Neutrinos) |l50l | is a "hybrid" (3(3 and solar experiment with ^°"Mo. It 
aims at studies of {mpp) with a sensitivity near 30 meV by measuring Oi/(3P decays of ^"^"Mo and the charged current 
^Be solar with an accuracy of about 10% by inverse () decays of ^^'^Mo. The decays to the ground and excited 
states are measured in prompt coincidence for the Oz//3/3 studies. The large Q=3.034 MeV results in a large phase- 
space factor to enhance the 0^(3(3 rate and a large energy to place the Ov(i(i signal above most background. MOON 
is a spectroscopic study of two /3 rays. As such, its capability to measure the energy and angular correlations for the 
two P rays can help identify the Oj//3/3 mechanism. It capability for a tight localization of (3(3 events in space and time 
is crucial for selecting Qv(3(3 and reducing background. 

A possible configuration of the MOON apparatus is a super module of hybrid plate and fiber scintillators with ^""Mo 
isotope mass totaling about 1 ton. One module of this apparatus consists of a thin (20 mg/cm^) Mo film interleaved 
between X-Y fiber planes and a plate scintillator. The fiber scintillators enable one to get the position resolution 
(1/K < 10~^ pixels per ton) and the scintillator plate provides an energy resolution (ct «2.2% at 3 MeV). A different 
detector option under consideration consists of foils within liquid Ar to obtain better energy resolution. Research 
has shown (1) that for the small plate scintillator, a = 2.2% , FWHM: 5% for Qv(3(3 including light attenuation, for 
MOON, (2) the position resolution of 3.2 x 10^, and (3) the feasibility of using centrifugal separation of MoFg gas 
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to produce -'^°°Mo enrichment of 85-90% in ton-scale quantities. Currently a proto-type MOON (MOON 1) is under 
construction (2003-2005). A proposal for MOON is planned by 2006. 

The MOON sensitivity can be evaluated as follows. The source is 1 ton of Mo with 85%^°"Mo. The 0i^f3f3 efficiency 
after energy and angle cuts is 0.28. The background from 21^/3/3 in the Oz//3/3 window during 5 years after cuts is 5.5(42) 
events for an energy resolution of (7=2.2(3)% (FWHM 5(7)%). Here a=2.2% is estimated based on the R&D data from 
a small-scale prototype and ELEGANT V, while 3% is a conservative value including a 50% factor. The background 
from cosmogenic isotopes and natural isotopes are less than 0.5 events. Then the 0vP(3 yield of ^/ background is 
2.3(6.4), corresponding to a half life 22 (7.7 )xlO^^. The mass sensitivities are, respectively, 13 (22) meV. Here the 
matrix element of |Afoi/| = 3 is estimated by referring to the recent 3 calculations. The MOON apparatus has the 
/3/3source separated from the detector, and therefore can be used for other /3/9 isotopes such as ^^Se, ^^°Nd and ^^^Cd 
as well by replacing Mo isotopes with other isotopes. 



Majorana 

The Majorana Collaboration proposes to field 500 kg of 86% enriched Ge detectors 14^. By using segmented 
crystals and pulse-shape analysis, multiple-site events can be identified and removed from the data stream. Internal 
backgrounds from cosmogenic radioactivities will be greatly reduced by these cuts and external 7-ray backgrounds 
will also be preferentially eliminated. Remaining will be single-site events like that due to f3p. The sensitivity is 
anticipated to be 4 x 10^^ y. 

Several research and development activities are currently proceeding. The collaboration is building a multiple-Ge 
detector array, referred to as MEGA, that will operate underground at the Waste Isolation Pilot Plant (WIPP) near 
Carlsbad, NM USA. This experiment will investigate the cryogenic cooling of many detectors sharing a cryostat in 
addition to permitting studies of detector-to-detector coincidence techniques for background and signal identification. 
A number of segmented crystals are also being studied to understand the impact of segmentation on background 
and signal. Recently, these segmented detectors have shown that the radial information provided by PSD is indeed 
orthogonal to the azimuthal information provided by the segmentation. The segmentation studies, or SEGA program, 
consists of one 12-segment enriched detector and a number of commercially available segmented detectors. Presently, 
commercially available segmented detectors are fabricated from n-type crystals. Such crystals are much more prone 
to surface damage and thus more difficult to handle when packaging inside their low-background cryostats. Hence 
the collaboration is also experimenting with segmenting p-type detectors. 

The Majorana design uses Ge detectors within a low-mass, electroformed Cu cryostat. Electroformed Cu is very 
free of radioactive contaminants. However, just how radio-pure the Cu is remains unknown. Hence, the collaboration 
plans to form Cu underground and study its radiopurity to a sensitivity below previous limits. 



Bare Ge Crystals 

The GENIUS collaboration^^ proposed to install 1 1 of enriched bare Ge crystals in liquid nitrogen. By eliminating 
much of the support material surrounding the crystals in previous experiments, this design is intended to reduce 
backgrounds of external origin. Note how this differs from the background-reduction philosophy ass ociat ed with 
pulse-shape analysis coupled with crystal segmentation. The primary advocates for this project indicate |l03l | th at its 
motivation has been questioned by their own claim of evidence for Oz^/3/3decay. Even so, the GENIUS test facilityll46| 
is being operated to demonstrate the effectiveness of operating crystals naked in liquid cryogen. 

Another group at the Max Plank In stitu te in Heidelberg, however, is proposing to pursue a similar idea. They have 
recently submitted a Letter of Intent flSlj to the Gran Sasso Laboratory. They propose to collect the enriched Ge 
crystals from both the Heidelberg-Moscow and IGEX expeiments and operate them in either liquid nitrogen or liquid 
argon. As a second phase of the proposal, they plan to purchase an additional 20-kg of enriched Ge detectors (most 
likely segmented) and operate with a total of 35 kg for about 3 years. Finally, they eventually plan to propose a large 
ton-scale experiment. It should be noted that this collaboration and the Majorana collaboration are cooperating on 
technical developments and if a future ton-scale experiment using ^^Ge proceeds these two groups will most likely 
merge and optimally combine the complementary technologies of bare-crystal operation and PSA-segmentation. 
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Nanocrystals 

Some elements may be suitable for loading liquid scintillator with metallic-oxide nanoparticles. Since Rayleigli 
scattering varies as the sixth power of the particle radius, it can be made relatively small for nanoparticles of radii 
below 5 nanometers. Particles of this size have been developed and commercial suppliers of Zr02, Nd2 03 etc. are 
available. Absorption of the materials must also be taken into account, but some of the metal oxides such as Zr02 
and Te02 are quite transparent in the optical region because of the substantial band gaps in these insulators. Some 
members of the SNO collaboration Tsij have been studying a configuration equivalent to filling the SNO cavity with a 
1% loaded liquid scintillator or approximately 10 t of isotope after the present heavy water experiment is completed. 
The group is currently researching the optical properties of potential nano-crystal solutions. In particular, one must 
demonstrate that sufficient energy resolution is achievable with liquid scintillator. 

Super-NEMO 

The currently operating NEMO-3 detector uses a tracking-calorimeter technique to detect /3/3. The source, is not 
the detector, but rather a thin foil located in the middle of a tracking chamber. This tracking chamber, which is made 
of geiger-drift cells, is surrounded by a calorimeter (scintillators with low radioactivity PMTs). The first of two long 
runs to search for Oi//3/3 decay began in February of 2003 with 7 kg of ^°'^Mo and 1kg of ^^Se. Additionally, 2 kg of 
various foils were placed in the detector to study 2i//3/3 decay and backgrounds. The first run will last approximately 
five years. So far a comprehensive study of the effects of various shields has been undertaken. The second run, again 
for five years, is currently planned to operate with 20 kg of *^Se. The *^Se is of particular interest because the 21^/3/3 
decay lifetime is 10 times that of ^°°Mo and thus may contribute less as a background to OvPP decay. 

There has been an unexpected background from radon, which is currently being fixed with a radon tent and a 
"radon trapping factory". Data collection free of the radon background will start in Fall 2004. The energy resolution 
is 250-keV FWHM at 3 MeV. Assuming the radon is eliminated with the radon tent and an efficiency of 14%, the 
background will be 5x10-* counts/keV/kg/yr. This yields T^J^ > 5 x lO^** in five years with (m^^) < 200-500 meV. 
During the second run with 20 kg of '^^Se, an efficiency of 14% and a background of 5x10^^^ counts/keV/kg/yr, the 
expected results in five years should be > 3 x 10^^ yr and (m^^) < 100 — 200 meV. 

The recent progress of the NEMO-3 Drogram (153i| has culminated in excellent 2vf3(3 results. In particular, the 
energy spectra from *°'^Mo contain approximately 10^ events and are nearly background free. These data permit, for 
the first time, a precise study of the spectra. In fact, there is hope that the data (Sutton 2004) will de mons trate 
whether the 2vl3l3 transition is primarily through a single intermediate state or through a number of states [l5^ . 

A much bigger project is currently being planned that would use 100 kg of source. The apparatus would have a 
large footprint however and the Frejus tunnel where NEMO-3 is housed would not be large enough to contain it. 
Currently the collaboration is studying the design of such a detector. 

XMASS 

The XMASS collaboration |l55[ Il5 9^ plans to build a 10 t natural Xe liquid scintillation detector. They expect an 
energy resolution of 3% at 1 MeV and hope to reach a value for T^j2 ^ ^■'^ ^ ^^'^'^ y- This detector would also be used 
for solar-neutrino studies and a search for dark matter. 

Borexino CTF 

In August of 2002, operations at the Borexino experiment resulted in the spill of scintillator. This led to the 
temporary closure of Hall C in the Gran Sasso Laboratory and a significant change in operations at the under- 
ground laborato ry. As a result, efforts to conv ert the Counting Test Facility (CTF) or Borexino itself into a Qvfip 
experiment , I37I Il53 | have been susDended jl6nl |. 
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How many Experiments are Required? 

In view of the importance and scale of new generation 0?^/3/3experiments, internationally cooperative efforts in both 
experiment and theory arc quite important. Thus it is reasonable to suggest a concentration of the limited resources 
on a few /3/3experiments However, it is critical to use different experimental techniques and isotopes to demonstrate 
the effect has really been seen and also to extract the most critical physics conclusions. It is therefore necessary to 
build a number of detectors and the reasons are enumerated here. 

1. The observation of a statistically significant signal in a single experiment might not be considered a discovery 
without clear confirmation from other independent experiments utilizing different isotopes. 

2. A nuclear matrix clement is necessary to deduce (m^^)from a measured Oi^/3/3rate. Since theoretical calculations 
of |Moi/|may include a substantial uncertainty, one needs experiments on different isotopes to extract a reliable 
value for the effective mass. 

3. Although light-neutrino exchange is the most natural explaination for Ovfipif it exists, there are other pos- 
sibilities. The relative matrix element values for different nuclei depend on the mechanism. Furthermore the 
matrix element situation is encouraging and one can anticipate a great improvement in the calculation precision. 
Therefore, measurements in several nuclei might be the most straight-forward way to provide insight into the 
mechanism of Oup(3. 

4. There are a number of different techniques being proposed for future experiments. Each has been previously used 

as effective prototypes for the proposals and therefore remains a strong candidate for future effort. However, 
the sensitivity of any given proposal depends strongly on its background estimate. It remains to be seen 
which of the technologies will successfully attain the required background. In addition, certain technologies 
provide capabilities such as measurments of the opening angle, individual electron energies, or the daughter 
production. These will not only help understand and remove background but they may also provide insight into 
the mechanism of Ov/Sp. 

VII. CONCLUSIONS 

Study of the neutrinoless double beta decay and searches for the manifestation of the neutrino mass in ordinary 
beta decay are the main sources of information about the absolute neutrino mass scale, and the only practical source 
of information about the charge conjugation properties of the neutrinos. Thus, these studies have a unique role in the 
plans for better understanding of the whole fast expanding field of neutrino physics. 

In this report we summarized the various aspects of the problem. We explained first the relation of the information 
that can be obtained from the analysis of and (3 decay experiments with the parameters describing neutrino 

oscillations and the absolute neutrino mass scale. We then discussed the nuclear structure issue and the uncertainty 
in determining the neutrino mass from a measured Oj//3/3 decay rate. We also briefly discussed the role of neutrino 
mass in cosmology, and the corresponding constraints on neutrino masses based on astrophysical observations. 

The remaining part of the report was devoted to the description of the existing and planned experiments. We 
described the present situation in Oi//3/3 and tritium beta decay, and in particular the recent, so far unverified, claim 
of the Oi^PP decay discovery. We concentrated then on the future plans for both areas of research. 

The situation in the direct neutrino mass searches using tritium f3 decay is relatively simple. There is only one 
realistic plan a very large new beta spectrometer which is being built in Germany. This KATRIN experiment has a 
design sensitivity approaching 200 meV. If the neutrino masses are quasi-degenerate, as would be the case if the recent 
double-beta decay claim proves true, KATRIN will see the effect. Although KATRIN is predominately a European 
effort, there is significant US participation. The design and construction of this experiment is proceeding well and we 
enthusiastically recommend the continuing strongly support of this program. 

There are many proposals, and even more ideas, for much larger Oi/(3p decay experiments than the existing ones. 
We outlined and justified in the report the strategy that we believe should be followed: 

1. A substantial number (preferably more than two) of 200-kg scale experiments (providing the capability to make 

a precision measurement at the quasi-degenerate mass scale) with large US participation should be supported as 
soon as possible. We estimate that the timescale of such experiments is 3-5 years, and that each such experiment 
will cost approximately $10M-$20M. 
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2. Concurrently, the development toward '--^l-ton experiments {i.e. sensitive to y^Am^^) should be supported, 
primarily as expansions of the 200-kg experiments. The corresponding plans for the procurement of the enriched 
isotopes, as well as for the development of a suitable underground facility, should be carried out. The US funding 
agencies should set up in a timely manner a mechanism to review and compare the various proposals for such 
experiments which span research supported by the High Energy and Nuclear Physics offices of DOE as well as 
by NSF. Each such experiment will cost approximately $50M-$100M and take 5-10 years to implement. 

3. A diverse R&D program developing additional techniques of Qv^fi decay study should be supported. 

The field of neutrino physics has made great strides recently. We believe that the study of 0v(3(i decay in particular 
could very well be the next one where significant discoveries will be made. 
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